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SUMMARY 


A fatigue-crack-propagation  model  is  developed  by  reducing  crack-tip 
plasticity  to  a single  pair  of  aggregate  dislocations,  termed  superdislocations, 
whose  response  is  equivalent  to  that  of  the  entire  plastic  zone.  The  crack-tip 
damage  or  deformation  and  accompanying  incremental  crack  advance  resulting  from 
each  loading  cycle  are  accumulated  in  a single  superdislocation  pair  for  resid- 
ual interaction  on  the  next  loading  cycle.  Key  features  of  this  development 
are  that  the  required  material  parameters  are  limited  to  the  shear  modulus,  G, 
Poisson's  ratio,  v,  and  tensile-yield  strength,  Y,  and  that  crack-surface 
closure  is  accommodated.  A significant  simulation  of  the  fatigue-crack-propaga- 
tion  process  along  with  a promising  correlation  with  experimental  results  are 
achieved. 

At  its  present  stage  of  development,  this  model  provides  a useful  charac- 
terization of  the  fatigue-crack-propagation  process.  Substantial  indications 
of  its  success  include 

(1)  A constant  crack-growth  rate  under  uniform  cyclic  loading. 

(2)  Predicted  crack-growth  rates  which  are  remarkably  close  to  the 
experimentally  determined  values  for  two  different  materials-- 
an  aluminum  alloy  and  a titanium  alloy.  The  predictions  for  the 
case  R = 0.5  are  particularly  outstanding. 

(3)  Predictions  from  fundamentally  correct  considerations  without 
the  use  of  empirical  parameters  of  any  kind.  The  only  parameters 
that  are  required  are  the  ordinary  mechanical  properties  G,  v, 
and  Y. 

(4)  An  automatic  inclusion  of  the  crack-closure  phenomenon. 

Areas  requiring  further  improvement  and  refinements  include  (1)  the  effect  of 
varying  the  stress  ratio,  R,  (2)  adequate  overload  perturbation  on  the  steady- 
state-growth  rate,  and  (3)  reduction  in  the  computational  time  required. 
Recommendations  for  the  next  stage  of  development  and  the  future  work  required 
are  presented  in  the  report. 
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1.  INTRODUCTION 


Fatigue-crack-propagation  rates  can  be  predicted  with  reasonable  accuracy 
using  empirical  relations  when  simple  constant-amplitude  cyclic  loading  is 
applied.  The  load  interaction  phenomenon  occurring  in  variable-amplitude  load- 
ing has  not  yet  been  properly  analyzed,  however.  This  precludes  the  accurate 
prediction  of  fatigue-crack-growth  rates  under  the  spectrum  (service)  loadings 
that  act  on  aircraft  and  other  engineering  structures.  To  provide  a quantita- 
tive formulation  for  fatigue-crack  growth  under  arbitrary  load  sequences,  a 
fundamentally  sound  analysis  technique  must  be  developed.  At  the  same  time, 
to  be  useful  for  routine  engineering  analyses,  the  required  analysis  procedure 
must  remain  as  simple  as  possible. 

The  goal  of  the  research  described  here  is  to  develop  an  accurate  computa- 
tional procedure  for  the  prediction  of  fatigue-crack  propagation  that  can  be 
used  in  engineering  applications.  To  achieve  this  goal 

(1)  The  model  must  be  capable  of  handling  load  cycles  that  vary 
arbitrarily  from  cycle  to  cycle  while  taking  the  load  inter- 
action effects  properly  into  account. 

(2)  The  material  properties  required  by  the  model  in  calculating 
crack-growth  rates  for  given  cyclic  load  histories  must  be 
based  upon  well-established  mechanical  properties  that  are 
independent  of  the  particular  load  spectrum  under  considera- 
tion. 

(3)  The  computational  procedure  evolved  must  be  efficient  enough 
to  enable  calculations  to  be  carried  out  over  load  histories 
corresponding  to  actual  service  conditions. 

To  meet  these  criteria,  a unique  approach  leading  to  a mathematical  predictive 
model  for  fatigue-crack  propagation  based  on  the  "inclined-strip-yield-super- 
dislocation"  representation  of  crack-tip  plasticity  has  been  devised.  The 
model  represents  an  effective  compromise  between  approaches  that  sacrifice 
some  of  the  basic  aspects  of  fatigue-crack  propagation  in  the  interest  of 
simplicity  and  a completely  rigorous  (and,  consequently,  exceedingly  cumber- 
some) approach. 

The  basis  of  the  model,  together  with  some  verifications  obtained  by 
comparison  with  established  results  in  fracture  and  fatigue,  is  given  in  this 


report.  Also,  an  experimental  program  on  Ti-6A1-4V  designed  to  implement  the 
development  of  the  model  and,  ultimately,  to  verify  its  predictions  is  de- 
scribed. Comparisons  between  the  predictions  of  the  model  and  experimental 
results  are  then  given.  Finally,  the  present  state  of  work  is  summarized  and 
recommendations  for  further  work  given.  For  the  convenience  of  the  reader, 
the  bulk  of  the  mathematical  and  experimental  details  have  been  omitted  from 
the  main  body  of  the  report  and  appear  in  appendices. 

2.  GENERAL  DESCRIPTION  OF  THE  MATHEMATICAL  MODEL 

There  is  currently  no  complete  agreement  on  the  fundamental  mechanisms 
that  control  fatigue-crack  propagation.  Nevertheless,  there  are  several 
features  of  the  problem  that  are  generally  accepted  as  playing  an  essential 
role  in  the  process.  These  must  be  taken  into  account  in  a proper  analysis. 
They  include 

• The  interaction  between  the  plastic  deformation  produced  by  the 
current  load  and  the  residual  plasticity  created  in  previous 
load  cycles. 

• The  connection  that  exists  between  the  local  deformation  at  the 
crack  tip--as  characterized  by  the  crack-tip-opening  displace- 
ment--and  the  crack-growth  increment. 

• The  closure  of  the  crack  faces  during  the  unloading  portion  of 
the  load  cycle  prior  to  reaching  the  minimum  load  level. 

While  analytical  models  have  been  developed  from  one  or  another  of  these  basic 
effects,  no  previously  developed  analysis  rigorously  incorporates  them  all. 
This  includes  models  employing  a finite-element  method  which  are  impractical 
beyond  a few  load  cycles. 

The  fatigue-crack-propagation  model  described  in  this  report  is  based  on 
the  well-established  fact  that  macroscopic  plasticity  can  be  considered  in 
terms  of  dislocation  arrays.  Figure  2.1  illustrates  how  this  basic  concept 
has  evolved  into  the  fatigue  model  being  developed  here.  The  following 
discussion  elaborates  on  the  conceptual  picture  given  by  Figure  2.1  and,  in 
addition,  delineates  the  individual  roles  played  in  this  work  by  the  various 
components  (e.g.,  the  superdislocation  concept)  of  the  complete  cycle-by- 
cycle fatigue-crack-propagation  model. 
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Plane  strain  plastic  deformation 
at  the  tip  of  a crack  under  fixed 
load 


Representation  of  crack-tip 
plasticity  by  dislocation  arrays 


Representation  of  crack-tip 
plasticity  by  a superdislocation 
pair  confined  to  slip  planes  ema1 
nating  from  the  crack  tip 


Representation  of  crack-tip 
plasticity  during  fatigue  by 
superdislocations 


Representation  of  crack-tip 
plasticity  during  fatigue  by 
combination  of  superdislocations 
and  super-superdislocations 


FIGURE  2.1 


EVOLUTION  OF  A FATIGUE -CRACK-GROWTH  MODEL  USING 
THE  INCLINED-STRIP-YIELD  SUPERDISLOCATION 
REPRESENTATION  OF  CRACK-TIP  PLASTICITY 
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Figure  2.1(a)  shows  a typical  plastic  enclave  surrounding  a crack  tip 
under  plane-strain  conditions.  Figure  2.1(b)  shows  an  equivalent  way  of 
characterizing  crack-tip  plasticity  through  the  use  of  dislocation  arrays.  It 
might  be  noted  that  generally  there  is  no  significant  advantage  of  the  disloca- 
tion point  of  view  over  the  conventional  continuum  approach.  However,  the 
dislocation  concept  can  be  reduced  to  a simpler  way  of  looking  at  the  problem 
that  can  offer  significant  computational  advantages  in  applications  to  fracture 
and  fatigue.  In  particular,  by  representing  crack-tip  plasticity  by  planar 
(one-dimensional)  dislocation  arrays  in  a "strip-yield-zone"  model,  consider- 
able mathematical  simplification  can  be  obtained  without  sacrificing  too  much 
accuracy.  This  fact  is  the  prime  motivation  for  adopting  a dislocation-based 
approach  here. 

Strip-yield-zone  characterizations  of  crack-tip  plasticity  are  most  easily 
made  if  the  dislocation  array  is  confined  to  the  plane  of  the  crack.  While 
reasonably  accurate  for  the  plane-stress  conditions  existing  in  a thin  section, 
this  is  a poor  representation  for  plane-strain  conditions  where,  as  indicated 
in  Figure  2.1(a),  plasticity  tends  to  spread  out  in  the  direction  normal  to  the 
crack  plane.  To  characterize  the  latter  case,  Bilby  and  Swinden  [1]*  took  the 
dislocation  slip  plane  to  be  inclined  to  the  crack  plane.  But,  they  were  then 
able  to  obtain  only  an  approximate  numerical  solution.  The  work  of  Atkinson 
[2-4]  circumvented  the  mathematical  difficulties  in  the  "inclined-strip-yield" 
model  by  introducing  the  idea  of  a superdislocation.  As  shown  in  Figure  2.1(c), 
the  superdislocation  is  considered  to  be  a dislocation  of  arbitrary  strength 
that  represents  the  net  effect  of  the  entire  dislocation  array  on  the  given 
slip  plane.  By  the  superdislocation  approach,  the  problem  is  reduced  to  two 
algebraic  equations  in  two  unknowns.  It  is  important  to  recognize  that, 
although  this  simple  representation  has  its  limitations  (e.g. , in  predicting 
the  plastic  zone  size),  it  offers  a very  accurate  prediction  of  the  crack-tip 
opening  displacement.  This  will  be  shown  later  in  this  report. 

The  mathematical  simplicity  of  the  superdislocation  representation  of 
crack-tip  plasticity  combined  with  an  accurate  prediction  of  the  state  of 
deformation  at  the  crack  tip  gives  this  approach  a strong  appeal  for  use  in  a 
fatigue-crack-propagation  model.  In  fact,  this  is  the  basis  of  the  model 
described  in  this  report.  That  is,  the  several  crack-tip-connected  phenomena 

* References  are  listed  on  page  65. 
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that  are  generally  accepted  as  playing  an  essential  role  in  the  fatigue-crack- 
growth  process,  namely  (1)  the  interaction  between  the  plastic  deformation 
produced  by  a load  and  the  residual  plasticity  created  in  previous  load  cycles; 
(2)  the  connection  that  exists  between  the  deformation  at  the  crack  tip  and  the 
crack-growth  increment;  and  (3)  the  fact  that  the  crack  faces  can  impinge 
during  the  unloading  portion  of  the  load  cycle  prior  to  reaching  the  minimum 
load  level  can  all  be  taken  into  account  with  the  superdislocation  on  an 
inclined-slip-plane  idea. 

At  this  point,  it  will  be  useful  to  identify  and  distinguish  between  the 
two  complementary  computational  approaches  that  are  being  pursued  in  this  work. 
The  first  is  the  complete  or  cycle-by-cycle  fatigue-crack-growth  model  in  which 
one  superdislocation  pair  is  generated  and  retained  in  the  computation  for  each 
load  cycle  considered.  This  is  shown  in  Figure  2.1(d).  The  second  approach 
is  the  super-superdislocation  idea  shown  in  Figure  2.1(e).  In  the  latter  case, 
the  residual  plasticity  is  represented  by  lumping  two  or  more  superdislocations 
together  to  reduce  the  computational  effort  required.  While  the  intention  is 
to  integrate  them  eventually,  they  are  being  pursued  separately  at  present. 
Therefore,  their  relative  advantages  and  disadvantages  will  be  discussed  here 
individually. 

The  complete  cycle-by-cycle  fatigue-cracx-growth  model  employing  the 
inclined  strip-yield-superdislocation  representation  of  crack-tip  plasticity 
has  now  been  substantially  developed.  This  model  is  believed  to  be  uniquely 
capable  of  providing  the  basis  for  meeting  the  objectives  of  this  research  as 
stated  in  the  Introduction.  Substantial  indications  that  the  model  will  be 
successful  already  exist.  As  described  in  more  detail  later  in  this  report, 
these  are  as  follows: 

(1)  After  a transient  period,  the  model  produces  a constant  crack- 
growth  rate  under  uniform  cyclic  loading.  Because  there  is 
nothing  built  into  the  model  that  will  force  this  to  happen, 
this  can  be  taken  as  one  verification  of  the  procedure. 

(2)  The  predicted  crack-growth  rates  for  the  uniform  cyclic  load- 

ing are  remarkably  close  to  the  experimentally  determined 
values  for  two  different  materials:  an  aluminum  alloy  and  a 

titanium  alloy.  The  predictions  for  the  case  R = 0.5  are 
particularly  outstanding. 
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(3)  The  model  can  be  used  to  make  predictions  from  fundamentally 
correct  considerations  without  the  use  of  empirical  param- 
eters or  "fudge  factors"  of  any  kind.  The  only  parameters 
that  are  required  are  the  ordinary  mechanical  properties  G, 
v,  and  Y. 

(4)  The  much-discussed  crack-closure  phenomenon  is  automatically 
taken  into  account  in  the  model.  This  is  done  both  simply  and 
accurately  without  requiring  any  ad  hoc  empirical  approxima- 
tions. This  is  accomplished  via  the  singularity  canceling 
equation  which,  when  satisfied,  assures  that  the  crack  forces 
will  close  smoothly  at  the  crack  tip. 

However,  the  present  cycle-by-cycle  model  suffers  from  the  following  disadvan- 
tages : 

(1)  The  effect  of  increasing  R while  holding  AK  constant  causes 
a modest  decrease  in  the  predicted  growth  rate.  The  experi- 
mentally observed  result  is  a marked  increase  in  growth  rate. 

(2)  While  an  overload  cycle  perturbs  the  steady-state  growth 
rate  for  a time,  the  observed  crack-growth  retardation 
phenomenon  in  which  the  crack-growth  rates  can  be  reduced 
by  two  orders  of  magnitude  is  not  reproduced. 

(3)  The  computer  time  required  to  perform  many  successive  load 
cycles,  while  not  unduly  costly  for  research  purposes,  is 
prohibitive  for  engineering  use. 

While  each  of  these  may  appear  to  be  a formidable  barrier  to  the  successful 
development  of  the  model,  a number  of  alternatives  have  been  identified  which 
may  provide  a means  of  circumventing  them.  These  are  described  in  the  conclud- 
ing section  of  this  report. 

The  use  of  the  inclined-strip-yield-superdislocation  representation  of 
crack-tip  plasticity  in  a mathematical  predictive  model  for  fatigue-crack 
propagation  has  suggested  the  use  of  a simpler  lumped  or  super-superdislocation 
model.  With  this  approach,  the  following  verifications  have  been  obtained. 

(1)  A prediction  of  the  crack-tip  crack-opening  displacement  under 
fixed  load  that  compares  very  well  with  highly  precise  finite- 
element  computations. 

(2)  A prediction  of  a steady-state  fatigue-crack-propagation  rate 
that  is  in  good  qualitative  agreement  with  experimental  results. 


(3)  A prediction  of  the  crack-growth  retardation  in  a uniform 
cyclic  loading  interrupted  by  a peak  overload  that  approxi- 
mates the  observed  effect  reasonably  well. 

In  the  first  of  these  verifications,  there  is  no  difference  between  the  super- 
dislocation  and  the  super-superdislocation  models.  In  the  latter  two  cases, 
the  super-superdislocation  represents  the  aggregate  effect  of  the  residual 
plasticity  associated  with  a multiplicity  of  superdislocations  that  were 
created  over  a large  number  of  load  cycles. 

It  may  be  useful  to  further  amplify  the  distinction  between  the  computa- 
tional models  based  on  the  superdislocation  and  super-superdislocation  con- 
cepts that  are  used  in  the  work  described  in  this  report.  The  ultimate  goal 
of  the  work  is  to  develop  a computational  model  that  can  be  applied  to  predict 
fatigue-crack-growth  rates  in  real  service  conditions.  This  calls  for  a 
mathematical  model  that  is  versatile  enough  to  accept  essentially  random 
changes  in  the  cycle-to-cycle  applied  loadings  while  remaining  simple  enough 
to  enable  it  to  be  applied  to  lengthy  load  histories.  The  work  reported  here 
has  progressed  far  enough  to  make  it  clear  that  a cycle-by-cycle  superdisloca- 
tion model  can  satisfy  the  first  of  these  requirements,  but  not  the  second. 
However,  a modification  to  lump  the  effect  of  possibly  ten  superdislocations 
into  one  or  more  single  degrees  of  freedom  super-superdislocations  can  allow 
the  model  to  achieve  the  second  criterion.  The  number  of  these  super-super- 
dislocations that  are  required  will,  of  course,  depend  upon  the  complexity  of 
the  load  spectrum  being  considered.  Nevertheless,  it  is  conceivable  that  the 
total  number  of  degrees  of  freedom  (superdislocations  plus  super-superdisloca- 
tions) need  not  exceed  ten.  If  borne  out  by  subsequent  work,  this  will  allow 
the  model  to  display  both  the  versatility  and  efficiency  that  it  must  have  to 
qualify  as  a basic  engineering  tool. 

While  the  prime  purpose  of  the  super-superdislocation  concept  is  as 
described,  it  has  also  been  put  to  a secondary  use.  The  secondary  use  is  in 
constructing  simple  models  for  steady-state  fatigue  with  and  without  interrup- 
tion by  an  overload.  These  models  are  not  intended  to  replace  the  cycle-by- 
cycle computational  model  in  any  sense.  Rather,  the  intent  is  to  have  a basis 
for  testing  the  various  concepts  and  computational  options  that  are  being 
considered  for  incorporation  into  the  complete  model.  The  linearized  super- 
superdislocation fatigue  models  permit  this;  first,  without  the  necessity  for 
using  computer  time  (because  simple  closed  form  relations  can  be  obtained  in 
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this  case)  and,  second,  because  an  interpretation  can  be  obtained  in  a clean- 
cut  manner.  Thus,  these  models  complement  the  complete  cycle-by-cycle  model. 

They  do  not  compete  with  it  but,  instead,  act  as  auxiliary  aids  to  assist  in 
its  development. 

To  help  make  this  essential  point  clearer,  the  above  description  is 
essentially  repeated  in  Table  2.1.  In  the  table,  each  level  of  the  model  is 
given  together  with  its  purpose,  its  strengths,  and  its  weaknesses  in  order  of 
ultimate  usefulness. 

3.  GENERAL  DESCRIPTION  AND  RESULTS  OF 
THE  EXPERIMENTAL  PROGRAM 

To  refine  and  confirm  a large  quantity  of  experimental  data  obtained  on 
previous  NADC-sponsored  programs  [5,6],  selected  fatigue-crack-propagation 
experiments  on  center-cracked  specimens  of  Ti-6A1-4V  alloy  plate  were  conducted 
to  provide  more  positive  substantiation  of  the  specific  analytical  developments 
of  this  research  program.  Included  were  the  generation  of  fatigue-crack- 
propagation  data  under  conditions  of  constant  cyclic  stress-intensity-factor 
range  (i.e.,  constant  AK) , the  study  of  the  overload  effects  on  crack-growth 
retardation,  the  calibration  of  specimen  compliance,  and  observations  of  crack- 
closure  effects. 

In  general,  for  basic  fatigue-crack-propagation  characterization,  data  are 
generated  under  constant-amplitude  loading.  However,  in  this  program,  to  permit 
a more  definitive  description  of  crack  behavior,  fatigue-crack-propagation  data 
were  generated  under  conditions  of  constant  cyclic  stress-intensity  factor  range, 
AK,  achieved  by  decrementing  the  load  with  increasing  crack  length  in  the  center- 
cracked  specimen.  This  approach  has  the  advantage  of  allowing  the  dependency  of 
fatigue-crack-propagation  rate,  da/dN,  on  K or  AK  to  be  more  clearly  established 
as  well  as  permitting  a large  number  of  overload  combinations  to  be  studied  over 
a range  of  crack  lengths  in  one  specimen. 

Both  mill-annealed  (£-inch  thick)  and  beta-processed  (5/16-inch  thick) 
Ti-6A1-4V  titanium  alloy  material  were  evaluated.  The  basic  experimental 
details  of  this  research  program  are  presented  in  Appendix  B.  Descriptions 
of  test  specimens,  materials,  and  experimental  details  of  this  research  pro- 
gram are  presented  in  Appendix  B.  Descriptions  of  test  specimens,  materials, 
and  experimental  procedures  are  included. 
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In  general,  the  results  were  quite  consistent  with  previously  generated 
data.  Of  surprise,  however,  was  the  lower  fatigue-crack-propagation-rate 
behavior  observed  in  the  beta-processed  material.  A comparative  summary  of 
these  and  previously  generated  rate  data  is  shown  in  Figure  3.1.  Other  sig- 
nificant results  included 

(1)  Fatigue-crack-growth  curves  generated  under  constant  stress- 
intensity-factor  range  (AK)  conditions  exhibited  a linear 
behavior  until  the  crack  exceeded  half  the  specimen  width 
confirming  the  dominance  of  AK  as  a primary  variable. 

(2)  The  beta-processed  material  exhibited  a much  coarser  macro- 
scopic surface  texture  than  the  mill-annealed  material  for 
equivalent  testing  conditions. 

(3)  Overloads  provided  definitive  retardation  (i.e.,  more  than 
10,000  cycles  of  retarded  growth)  only  when  the  overload 
exceeded  the  nominal  maximum  cyclic  stress  by  a factor  of 


2 (i.e.,  100  percent  overload). 

(4)  Compliance  response  of  the  specimen  was  consistent  with  the 
analytically  derived  values. 

(5)  Crack  closure,  although  elusive  to  evaluate  quantitatively, 
was  a very  detectable  phenomenon. 


4.  BASIC  DESCRIPTION  OF  THE  COMPUTATIONAL 
MODEL  AND  THE  ANALYSIS  PROCEDURE 


Preliminary  Discussion 


The  mathematical  derivation  of  the  basic  equations  for  the  cycle-by-cycle 
predictive  model  for  fatigue-crack  growth  under  arbitrary  spectrum  loading  is 
given  in  Appendix  A of  this  report.  For  this  reason,  it  will  only  be  out- 
lined here. 

The  essential  idea  upon  which  the  model  is  based  is  that  microscopic 
crack-tip  plasticity  generally  can  be  represented  by  dislocation  arrays.  To 
simplify  the  mathematics  for  fatigue-crack  propagation,  it  is  further  assumed 
that  the  plastic  deformation  created  by  a cycle  of  load  can  be  adequately 
represented  by  a single  superdislocation.  That  is,  the  plastic  deformation 
symmetric  with  respect  to  the  crack  plane  is  represented  by  a single  pair  of 


Previous  Program 


Test  data  bands 


Rate  computed  by 
Equation  B-4 


Present  Program 


Material 


Beta  Mill  Stress 
Processed  Annealed  Ratio,  R 


Stress  Intensity  Range,  AK,  ksi  — in 


FIGURE  3.1.  COMPARISON  OF  PRESENT  AND  PREVIOUS  DATA 
(See  Appendix  B also) 


superdislocations  symmetrically  placed  with  respect  to  the  crack  plane.  A 
final  major  assumption  gives  rise  to  the  super-superdislocation  model;  namely, 
that  the  residual  plasticity  created  in  a series  of  load  cycles  can  be  aggre- 
gated into  one  or  more  super-superdislocations.  This  process  has  been  illus- 
trated in  Figure  2.1. 

Modifications  to  the  mathematical  model  described  in  Section  2 have  been 
made  and  can  be  mentioned  here.  The  first  modification  was  instituted  because 
of  a difficulty  in  the  cycle-by-cycle  model  when  the  minimum  load  was  set  to 
zero.  In  this  case,  all  residual  plasticity  could  vanish.  To  overcome  this 
objectionable  feature,  a second  pair  of  superdislocations  exists  on  the  crack- 
tip  slip  plane  on  the  minimum  load  portion  of  the  load  cycle.  This  has  the 
effect  of  locking  the  nascent  dislocations  in  place  when  the  applied  stress 
vanishes.  The  physical  basis  for  the  existence  of  the  second  pair  is  crack 
closure  which  represents  the  plasticity  due  to  the  contact  stresses  acting  on 
the  crack  faces  where  they  impinge. 

The  second  and  more  fundamentally  important  modification  is  connected 
with  the  role  of  the  singularity-canceling  equation  and  its  relation  to  the 
crack-closure  phenomenon.  In  the  beginning  stages  of  this  work,  the  notion 
that  the  crack-tip  singularity  must  be  abolished  (as  in  the  Dugdale  model)  was 
not  questioned.  Eventually,  two  inconsistencies  were  recognized.  First,  the 
dislocations  introduced  to  represent  plastic  deformation  have  singular  points 
that  remain  when  the  crack-tip  singularity  is  removed.  Second,  removing  the 
crack-tip  singularity  is  tantamount  to  setting  the  linear-elastic  fracture- 
mechanics  strain  energy-release  rate  parameter,  to  zero.  As  will  be  shown 
later  in  this  report,  Jt  has  the  physical  interpretation  of  a driving  force  for 
crack  advance,  allowing  crack  growth  with  Jf  equal  to  zero  is  inconsistent  with 
the  basic  notions  of  fracture  mechanics.  The  proper  role  for  the  singularity- 
canceling equation,  it  is  now  clear,  is  in  connection  with  crack  closure  under 
a reduced  load  in  a cyclic  loading  sequence. 

An  illustration  of  the  superdislocation  model  at  the  point  at  which  the 
first  minimum  load  is  imposed  in  a loading  sequence,  during  which  the  crack 
has  grown  a measurable  amount,  is  shown  in  Figure  4.1.  It  is  supposed  that 
the  load  sequence  has  been  such  that  crack  closure  has  occurred  necessitating 
the  use  of  the  singularity-canceling  equation.  The  crack-face  displacements 
would  then  be  roughly  as  shown  (NB,  the  figure  is  not  drawn  to  scale).  Please 
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applied  load 


note  that  the  plastic  zones  portrayed  in  this  and  the  following  figures  are 
not  computed  in  the  model.  They  are  shown  only  to  clarify  the  way  in  which 
the  superdislocation  idea  can  represent  crack-tip  plasticity. 

The  situation  when  the  load  is  increased  to  a maximum  level  is  illus- 
trated in  Figure  4.2.  Here  the  crack  tip  is  assumed  to  be  open  and  to  have 
the  characteristic  blunted  appearance  of  a crack  tip  surrounded  by  a plasti- 
cally deformed  region.  It  should  be  emphasized  that  the  open  appearance 
occurs  only  when  the  crack-tip  singularity  is  not  abolished  at  the  maximum 
load  level.*  Having  the  equilibrium  solution  for  this  state,  the  crack  advance 
increment  can  be  calculated,  the  superdislocation  adjusted  accordingly,  and 
the  load  reduced.  The  crack  tip  would  again  appear  as  in  Figure  4.1. 

Basic  Equations  of  the  Model 


The  approach  to  be  followed  in  developing  the  basic  equations  for  the 
model  is  completely  within  the  linear  theory  of  elasticity.  The  reason  is 
that  crack-tip  plasticity  is  represented  by  discrete  singularities  (disloca- 
tions). The  most  convenient  representation  is  the  complex  variable  method  of 
Muskhelishvili  in  which  the  stress  and  displacement  fields  are  given  by 

axx  + ayy  = 4 Re  t®7  (z)  1 

2 (CTyy  " axx)  + iTxy  = z 0"  (z)  + (z)  ’ (4,2) 

2G  (u  + iv)  = k<S(z)  - z 2 ' (z)  - i|r  (z)  . (4.3) 

where  G is  the  shear  modulus,  z is  the  complex  variable  given  by 

z = x + iy  , 

3-4v  , plane  strain 
3-v 

, plane  stress 

and  2 and  i|r  are  functions  of  a complex  variable  that  depend  on  the  boundary 
conditions . 


* Circumstances  could  be  envisioned  in  which  the  crack  would  remain  closed 
even  at  maximum  load.  It  is  possible  that  crack  arrest  following  a large 
overload  occurs  in  this  way. 
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FIGURE  4.2.  ILLUSTRATION  OF  SUPERDISLOCATION  REPRESENTATION  OF 
CRACK- TIP  PLASTICITY  SHOWING  CRACK  OPENING  AT 
MAXIMUM  LOAD  LEVEL  IN  A CYCLIC  LOADING  SEQUENCE 


The  technique  used  for  determining  the  solution  to  the  problem  of  interest 
employs  linear  superposition  of  potential  junctions  valid  for  an  infinite 
domain  that  satisfy  the  boundary  conditions  for  four  subproblems.  These  are 

(1)  Remotely  applied  uniform  stress  (e.g.,  uniaxial  tension, 
biaxial  tension,  shear)  for  infinite  domain  without  cracks 

(2)  A crack  with  surface  tractions  equal  and  opposite  to  the 
applied  stress 

(3)  Isolated  dislocation  pairs 

(4)  Crack  with  surface  tractions  equal  and  opposite  to  the 
stresses  arising  from  the  dislocation  pairs. 

Note  that  Subproblems  (2)  and  (4)  give  rise  to  singular  behavior  at  the  crack 
tip.  Because  the  form  of  the  singularity  is  the  same  in  each  case  (i.e.,  r“^), 
these  can  be  made  to  cancel,  if  desired. 

The  key  element  in  the  approach  is  the  use  of  a single  superdislocation 
pair  (in  the  case  where  the  deformation  is  symmetric  with  respect  to  the  crack 
plane)  to  represent  the  plastic  zone  in  each  load  cycle.  This  requires  only 
two  new  unknowns,  strength  and  the  position  of  a superdislocation,  to  be 
determined  in  each  load  cycle.  The  unknowns  are  determined  by  satisfying  a 
force  equilibrium  relation  for  a superdislocation  and  imposing  a singularity- 
canceling equation  at  the  crack  tip  to  insure  that  the  crack  closes  smoothly 
without  interpenetration  at  the  tip. 

Consider  that  a tensile  loading  is  applied  in  the  direction  normal  to  the 
crack  plane  and  that  plane-strain  conditions  exist.  Then  the  force  equilibrium 
equation  is 
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t = eh  + -r-r. r Zb.  (g.  +k.  ),n  = l,2,...M,  (4.5) 

n n 4n ( 1 - v ) j = 1 j °jn  jn 

while  the  singularity-canceling  or  crack-closure  equation  is 


OG  — / . ■»  , b . f . , 

4rr(l  - v)  j j 

where  a = applied  tensile  stress 
a = half  crack  length 
M = number  of  load  cycles  imposed 

t h 

bj  = strength  of  the  dislocation  created  in  the  j — load  cycle 


(4.6) 


22 


G = shear  modulus 
v = Poisson's  ratio  . 

with  the  remaining  undefined  quantities  in  Equations  (4.5)  and  (4.6)  being 
known  functions  of  the  dislocation  position  in  terms  of  the  complex  variable 
representation  given  by 


Z . = a . + i.e 
J J J 


j =1,2...  M , 


where  9 = angle  between  the  crack  plane  and  slip  plane  (constant) 

t h 

a . = crack  length  at  the  time  that  the  j dislocation  was  emitted 

^ th 

l = distance  from  the  crack  plane  along  the  slip  plane  to  the  j 

dislocation. 

Before  writing  out  these  expressions  explicitly,  it  may  be  useful  to  point  out 
their  physical  interpretation.  Note  that  the  symbol  =>  is  taken  to  mean  "is 
associated  with"  in  these  definitions 


f =>  contribution  to  the  singularity-canceling  equation  due  to  disloca- 
n 

tion  pair  generated  in  the  n load  application  (dimensionless). 

g.  =>  stress  acting  along  the  slip  plane  at  the  position  occupied  by  the 

th 

dislocation  generated  in  the  n load  application  (first  quadrant) 

due  to  the  surface  tractions  that  abolish  the  crack  face  stresses 

t h 

caused  by  the  dislocations  generated  in  the  j load  application 
(dimensions  of  reciprocal  length). 

h^  =>  stress  acting  along  the  slip  plane  at  the  position  occupied  by  the 

dislocation  generated  in  the  nC^  load  application  (first  quadrant) 

due  to  the  applied  stress  (dimensionless). 

k.  =>  stress  acting  along  the  slip  plane  at  the  position  occupied  by  the 

th 

dislocation  generated  in  the  n load  application  (first  quadrant) 
due  to  the  stress  field  of  the  dislocations  generated  in  the  jCh 
load  application  (dimensions  of  reciprocal  length). 

Omitting  the  details  involved  in  their  derivation,  the  quantities  defined 
on  the  preceding  page  are  found  to  be  as  follows: 


f = 4 sin  9 Re 
n 


+ 2 i Re  7 v,  i 

° \«n-^Zn  - = / 


h = sin  0 cos  9 + or  l Re 
n n 


(4.7) 


(4.8) 
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The  manner  in  which  individual  superdis location  pairs  are  assembled  into  a 
complete  cvcle-by-cycle  model  governed  by  the  above  equations  is  illustrated 
in  Figure  4.3. 

It  is  considered  that  the  material  provides  an  intrinsic  resistance  to 
the  movement  of  the  superdislocation.  This  parameter  is  denoted  r (for  inter- 
nal) and  is  related  to  Y,  the  tensile  yield  stress  by 

t t = v/s  . (4.ii) 

The  equations  of  force  equilibrium  that  are  the  basis  for  the  mathematical 
model  are  formed  by  setting  the  right-hand  side  of  Equation  (4.5)  equal  to  t 
for  each  of  the  M superdislocations  in  the  problem. 

An  expression  for  the  crack-opening  displacements  can  also  be  determined. 
Again  omitting  the  details,  the  result  is  found  to  be 


Av 


M 

sin  0 Z 

J-l 


(4.12) 


where  Av  is  the  crack-opening  displacement  at  a position  x <.2.  In  Equation 
(4.12),  x is  measured  from  the  center  of  the  crack  and  d^  is  the  function  given 
by 


Linearized  Model  for  Fixed-Load  Problem 


In  this  section,  the  problem  illustrated  by  Figure  4.4(a)  will  be  consid- 
ered. Because  a single  superdislocation  pair  suffices,  the  subscripts  will 
be  omitted.  The  notation  is  then  as  shown  in  Figure  4.4(b). 

The  governing  equations  are  for  the  force  equilibrium  of  one  of  the  dis- 
locations (only  one  needs  to  be  considered  because  of  symmetry)  and  the  singu- 
larity-calceling  equation.  The  unknowns  are  the  superdislocation  strength,  b, 
and  the  position,  l.  Setting  t = in  Equation  (4.5),  it  is  found  that 

= sh  (A)  + ~ (^VL)  CS  U>  + k («]  , (4.14) 


and 


a 


Gb 

TT  (<  + 1)  a 


f U)  . 


(4.15) 


In  the  case  of  small-scale  yielding  where  A « a,  Equations  (4.7)  - (4.10) 
reduce  to 


f(£)  = 6 sin  6 cos  j 

gU)  =■  - j (1  - cos  2 6) 

h(i)  = ~ 2 sin  ® cos  2 + sin  ® cos  ® 

k (A)  = - — cos  29 

AJ 


(4.16) 


Further  assuming  that  a « Y and  substituting  Equation  (4.16)  into  Equations 
(4.14)  and  (4.15)  gives  the  following  relation  for  the  superdislocation 
position: 

3 

l 1 L . „ e 1 l /a  \s 

a = 24  3 sln9cos  2 " — : r (y)  • 

< sin  9 cos  — ’ 

In  this  work,  the  particular  angle  5 selected  will  be  that  which  maximizes 
the  extent  of  yielding  (equivalently,  the  plane  of  maximum  shear.  From  Equation 
(4.17)  this  is 

27 


— 


X-  2 


Plastic  region 


Inclined  Strip*Yield-Zone-Superdislocation  Model  for  Small 
Scale  Yielding  (Symmetric  Tensile  Loading) 


-Slip  plane 
Superdislocation 


Crack  face 


(b)  Detail  of  Crack-Tip  Deformation 


FIGURE  4.4 


APPLICATION  OF  SUPERDISLOCATION  MODEL  TO  SMALL-SCALE 
YIELDING  UNDER  CONSTANT-APPLIED  STRESS 


which  turns  out  to  be 


de  i sin  9 cos  ! } m 0 


(4.18) 


0 = cos"1  (1/3)  = 70.53°  . 


(4.19) 


This  is  the  value  used  in  this  work. 

The  superdislocation  cannot  be  expected  to  give  precise  information  about 
the  size  and  shape  of  the  plastic  zone  per  se.  However,  this  will  be  unim- 
portant if  it  does  reflect  with  reasonable  accuracy  the  effect  of  the  plastic 
deformation  at  the  crack  tip.  Consider  the  prediction  of  6,  the  crack-tip- 
opening displacement  in  plane  strain  under  fixed  load  for  small-scale  yield- 
ing conditions.  This  can  be  written  as 

6 = 2b  sin  0 . 

Upon  solving  for  b and  substituting  into  the  above,  it  is  found  that 


6 = 


2 

3/5 


cos 


0/2 


3 sin  0 


cos  — - 


sin  0 cos  "r 


t- 

EY 


(4.20) 


where  K = cj/tt2 


Hence , 


(4.21) 


where  a ■ 0.433  for  0 = 70.53°. 

A comparison  between  the  results  obtained  from  a variety  of  both  approximate 
and  highly  rigorous  calculations  of  the  crack-opening  displacement  are  given 
in  Table  4.1.  It  can  be  seen  that  the  simple  superdislocation  model  compares 
very  favorably  with  the  more  precise  results. 


Relation  Between  the  Singularity-Canceling 
Equation  and  Crack  Closure 


Using  the  same  mathematical  procedures  that  led  to  the  singularity- 
canceling equation,  an  expression  can  be  devised  for  the  combined  strength  of 
the  singular  term.  Because,  in  general,  the  relation  between  the  strain- 
energy-release  rate,  and  K,  the  stress-intensity  factor  is  given  for  plane 
strain  by 
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TABLE  4.1.  RELATION  BETWEEN  STRESS  INTI’S ITY- FACTOR  AND  CRACK* OPENING 
DISPLACEMENT  FOR  SMALL-SCALE  YIELDING  CONDITIONS 


, K2 

6 = EY 

Investigators 

Technique 

a 

Wells 

Linear  Elasticity 

1.159 

Goodier  and  Field 

Dugdale  Model 

0.910 

Rice 

Deformation  Plasticity 
(J  Integral) 

0.613 

Rice  and  Johnson 

Deformation  Plasticity 
(J  Integral) 

0.717 

Levy,  Marcal,  Ostegren, 
and  Rice 

Incremental  Plasticity 
(Finite-element  method) 

0.425 

Rice 

Bilby-Swinden  Model 
(Approximate) 

0.440 

Rice  and  Tracey 

Incremental  Plasticity 
(Singular  Finite  Element) 

0.493 

This  report 


Inclined  Strip  Yield 
(Superdislocation  Model) 


0.433 


K6 


T \ 


l-v3  ’ 
it  is  found  that 


. 1 -v3 

J = — - — nj 


a - 


TT  (K  + 1)  a 


j-1 


b.  f. 
J J 


(4.22) 


Comparison  of  Equation  (4.22)  with  Equation  (4.12)  shows  that  the  bracketed 
term  is  just  equal  to  the  combined  coefficient  of  the  singular  term  and,  hence, 
is  set  equal  to  zero  to  abolish  the  crack-tip  singularity. 

In  fracture  mechanics,  Jt  has  the  physical  interpretation  of  being  the 
driving  force  for  crack  advance.  Consequently,  allowing  crack  growth  to 
occur  with  equal  to  zero  (as  would  be  the  case  if  Equation  (4.12)  is  imposed 
during  the  maximum  load  in  fatigue)  is  inconsistent  with  the  basic  notions  of 
fracture  mechanics.  A more  appropriate  role  for  the  singularity-canceling 
equation  may  be  in  connection  with  crack  closure  under  a nonzero  load  in  a 
cyclic  loading  -'equence.  This  can  be  seen  in  connection  with  the  crack  face 
displacements,  Equation  (4.10). 

Consider  small-scale  yielding  conditions  and  confine  attention  to  a 
point  on  the  crack  face  in  the  near  neighborhood  of  the  crack  tip.  Let  the 
distance  from  the  crack  tip  be  given  by  r = a - x.  Then,  introducing  the  condi- 
tions r « CL  and  i « a in  Equation  (4.10),  it  is  found  that 


M 


...  k ) <j(*+1)  'a\H  2 e „ . _ 

Av(r)  = r | — ^ — *-  ~ sin  6 cos  - E b^D^ 


(4.23) 


where  D 


V* 


2(1  - r) 


j i ,z  + 2ri . cos  9 + r3  l } - 2x1 . cos  9 + r3 
J J J J 

Now,  if  the  right-hand  side  of  Equation  (4.23)  is  expanded  in  powers  of  r, 
the  result  can  be  shown  to  be 

Av  (r)  = r*  j 2l£±Ii  (ijs  . 1 sin  9 cos  | J + 9 (r*)  . (4.24) 

Under  small-scale  yielding  conditions,  the  coefficient  of  the  r^  term  is  just 
the  singularity-canceling  equation*.  Consequently,  by  removing  the  crack-tip 


* To  see  this,  substitute  the  first  of  Equations  (4.16)  into  Equation  (4.6), 
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singularity,  the  crack-face  displacements  will  be  of  the  order  r2.  This  means 
that  the  slope  of  the  displacements  will  be  of  order  r2.  At  the  crack  tip, 
r * 0.  Hence,  satisfying  the  singularity-canceling  equation  forces  the  crack 
faces  to  close  smoothly  at  the  crack  tip  with  no  interpenetration. 

This  is  obviously  a necessary  condition  for  crack  closure  (although  not 
sufficient)  and,  therefore,  it  is  appropriate  to  impose  it  at  the  minimum  load 
in  a cyclic  loading  sequence.  Conversely,  it  is  inappropriate  to  impose  the 
singularity-canceling  equation  at  the  maximum  load  as  has  been  done  routinely 
in  all  previous  work.  To  do  so  is  inconsistent  with  all  observations  on 
fatigue-crack  growth  (which  reveals  that  the  crack  blunts  on  the  increasing 
load  part  of  the  cycle)  and  with  the  idea  that  crack  growth  can  take  place 
with  zero  J>.  This  condition  should  be  imposed  only  at  the  minimum  load,  there- 
fore. 

The  Crack-Growth  Criterion  and 
the  Solution  Procedure 

As  already  noted,  the  mathematical  simplicity  of  the  superdislocation 
representation  of  crack-tip  plasticity  combined  with  an  accurate  prediction  of 
the  state  of  deformation  at  the  crack  tip  gives  this  approach  a strong  appeal 
for  use  in  a fatigue-crack-propagation  model.  In  addition,  the  displacement 
field  associated  with  the  superdislocations  distorts  the  crack-opening  displace- 
ment (COD)  in  the  proximity  of  the  crack  tip,  allowing  the  possibility  of  crack 
closure  under  positive  load.  At  the  same  time,  an  intrusion  into  the  unbroken 
material  at  the  crack  tip  is  predicted.  This  can  be  taken  as  a measure  of  the 
crack  advance  increment  in  the  load  cycle.  The  latter  effect  arises  from  the 
displacement  discontinuity  (which  is  associated  with  any  dislocation)  inter- 
secting the  crack  at  its  tip,  giving  a "sliding  off"  region  there. 

Because  the  magnitude  of  the  displacement  discontinuity  is  just  equal  to 
the  superdislocation  strength,  the  crack-growth  rate  is  simply  related  to  the 
strength  of  the  most  recently  emitted  superdislocation.  Hence,  the  crack 
growth  at  the  NC^  cycle  of  loading  is 


(«)„  ■ bN  C°S  9 • 


(4.25) 


where  9 is  the  angle  between  the  superdislocation  slip  plane  and  the  crack 
plane.  Taking  the  slip  plane  to  be  the  plane  of  maximum  shear  gives  9 = cos-i 
cos-1 (1/3)  = 70.3°  as  determined  in  Equation  (4.19). 
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The  solution  procedure  involves  an  iterative  technique  to  determine  the 
superdislocation  strengths  and  positions  in  each  load  cycle.  This  is  done  by 
solving  the  system  of  nonlinear  algebraic  equations  given  by  Equations  (4.5) 
and  (4.6).  A Runge-Kutta  method  is  currently  used.  Having  the  solution,  the 
crack-growth  increment  is  obtained  from  Equation  (4.25).  It  can,  therefore, 
be  recognized  that  no  arbitrary  disposable  parameters  are  introduced  in  order 
to  determine  the  crack-growth  rates.  The  only  parameters  that  enter  the 
computation  are  the  ordinary  mechanical  properties  of  the  material  being 
considered  and,  of  course,  the  cyclic  load  history. 

The  iterative  procedure  used  to  determine  the  equilibrium  dislocation 
positions  considers  that  the  internal  fraction  stress,  acts  along  the  slip 

plane  to  always  oppose  dislocation  motion.  Therefore,  the  net  stress,  F , 

ch.  ^ 

acting  to  move  The  n dislocation  along  the  slip  plane  is  (NB,  positive  force 

acts  away  from  the  crack  plane,  negative  force  acts  towards  the  crack  plane) 


l rn  - sgn(Tn)Ti  , |rn|  s t. 

where,  as  above,  is  taken  to  be  a material  constant  related  to  the  yield 

stress  (i.e.,  t = Y/3)  and  sgn  (x)  = x/  |x|  . Note  that  the  friction  force 

resisting  movement  back  towards  the  crack  tip  need  not  (and  probably  is  not) 
equal  to  that  resisting  motion  outwards. 

In  detail,  the  procedure  used  in  determining  fatigue-crack-growth  rates 
under  arbitrary  cycle-to-cycle  loading  starting  from  an  undeformed  (no  previous 
loads)  crack  tip  is  described  next.  It  is  understood  that  the  material's 
elastic  modulus,  Poisson's  ratio,  and  yield  stress,  together  with  an  initial 
crack  length,  have  all  been  specified  in  addition  to  the  loading  history  to 
be  followed.  For  definiteness,  it  will  be  considered  that  the  latter  quan- 
tities have  been  given  in  terms  of  K and  K . values  for  each  load  cycle. 

max  mm  1 

Then,  the  computational  procedure  is  as  follows  (subscripts  denote  cycle 
numbers ) : 

(1)  Determine  the  position  along  the  slip  plane  occupied  by  the 

first  dislocation,  , from  Equation  (4.17)  using  a = (K  )j / 

(nag)2.  (Note  that  although  this  is  an  approximate  linearized 
version  of  the  complete  equation,  it  will  give  an  accurate 
enough  result  so  that  no  further  iteration  is  necessary.) 
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(2)  Keeping  the  applied  stress  the  same,  compute  the  superdisloca- 

tion strength,  bx , from  a linearized  version  of  the  singularity- 
canceling equation.  Equation  (4.6).  Then,  from  Equation  (4.25), 
calculate  the  crack-growth  increment  (A2)x  from  the  relation 
(A2)1  = bx  cos  6. 

(3)  Let  the  crack  extend  such  that  its  half  length  2 is 
2 = 2q  *+■  (A  2 )x 

where  20  is  the  initial  crack  length. 

(4)  Reduce  the  applied  load  to  cr  = (Kmin)i / 3) 2 and  allow  a second 
superdislocation  of  strength  bx  ' < 0 to  appear  on  the  slip  plane. 
Find  an  equilibrium  solution  with  the  crack-closure  condition 
satisfied.  Then  assume  that  as  the  load  is  once  again  increased, 
the  inner  dislocation  combines  with  the  outer  to  give  a super- 
dislocation of  strength  bx  + bx  ' at  the  position  L 1 . Note  that 
this  strength  remains  fixed  for  the  remainder  of  the  computation. 

From  this  point  on,  the  computational  procedure  is  the  same  regardless  of  the 
number  of  superdislocations  and  load  cycles  performed.  So,  consider  that  M - 
superdislocation  pairs  have  been  generated  in  M - 1 load  cycles  where  M > 1. 
The  computation  for  the  Mtb  load  cycle  then  proceeds  as  follows. 


(5)  Increase  the  applied  load  to  a = (K  )w/(tt2) 
rr  max  M 


and  consider  that 
a nascent  dislocation  exists  at  the  crack  tip.  Assume  that  all 
dislocations  move  along  the  slip  plane  at  a speed  that  is  propor- 
tional to  the  net  force  acting  on  them;  i.e., 
di. 

= F , j = 1,  2,  ...  , M . 

Use  the  Runge-Kutta  method  to  allow  the  superdislocation  to  move  to- 
wards their  equilibrium  positions.  The  strength  of  the  nascent 
superdislocation,  b^,  will  be  determined  by  satisfying  the 
singularity-canceling  equation. 

Continue  the  iterative  solution  procedure  until  equilibrium  has 
been  achieved.  Then,  calculate  the  crack-growth  increment 

(A2)M  = bM  C°S  9 ’ 
and  increase  the  crack  length  to 
M 

2 = 20  + I (A  2 ) . . 

j=l  J 


(6) 
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(7)  Reduce  the  applied  load  to  a = C17 a ) * allow  the  secondary 

superdislocation  b1  < 0 to  appear,  determine  an  equilibrium  solu- 
tion that  satisfies  the  crack-closure  condition,  then  replace  b 
by  b^  + b^1  at  This  completes  the  load  cycle. 

Steps  5 through  7 are  repeated  for  each  load  cycle  after  the  first. 

Table  4.2  shows  the  results  taken  from  a typical  computation.  In  this 

h * 

example,  a uniform  load  history  given  by  Kmax  = 18  ksi  in  , K = 9 ksi  in 
was  applied  to  a plate  with  30  = 1 inch.  The  material  being  simulated  is 
aluminum. 

Super-Superdislocation  Model  for 
Steady-State  Fatigue-Crack  Growth 
Under  Uniform  Cyclic  Loading 

In  a development  similar  to  that  given  in  a preceding  section,  a linear- 
ized version  of  the  model  will  be  used  to  consider  the  case  of  steady-state 
fatigue-crack  growth  under  uniform  cyclic  loading.  Specifically,  it  will  be 
considered  here  that  the  applied  loads  will  be  adjusted  for  the  increase  in 
crack  length  such  that  K and  K . will  be  the  same  for  each  load  cycle.  As 
in  the  above,  the  "damage"  accumulated  in  the  form  of  plastic  deformation 
during  a period  of  loading  will  be  represented  by  a single  lumped  dislocation. 
The  differences  are  now  that  (1)  it  will  be  assumed  that  the  lumped  dislocation 
moves  with  the  crack  tip  in  a steady-state  fashion  and  (2)  it  represents  the 
accumulation  over  a large  number  of  load  cycles  rather  than  just  one.  Conse- 
quently, while  it  was  appropriate  to  refer  to  the  lumped  dislocation  in  the 
single  fixed-load  situation  as  a superdislocation,  in  the  present  problem,  it 
must  be  called  a super-superdislocation. 

The  primary  virtue  of  the  super-superdislocation  approach  is  that  the 
number  of  variables  in  the  problem  is  considerably  reduced.  Very  economical 
solutions  are  therefore  possible.  But,  there  is  also  a price  to  be  paid  to 
achieve  this  simplicity.  One  such  disadvantageous  feature  is  that  the  close 
connection  between  the  crack-tip-opening  displacement  (CTOD)  and  the  strength 
of  the  superdislocation  emanating  from  the  crack  tip  is  lost.  Hence,  the  crack- 
growth  criterion  based  directly  on  the  CTOD  cannot  now  be  used.  Instead,  it 
will  be  assumed  that  the  crack-advance  increment  is  proportional  to  the  strength 
of  the  super-superdislocation  with  the  constant  of  proportionality  being 
inferred  from  results  of  the  cycle-by-cycle  computation. 
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TABLE  4.2.  EXAMPLE  COMPUTATIONAL  RESULT  FOR  FATIGUE-CRACK 
GROWTH  UNDER  UNIFORM  CYCLIC  LOADING  IN  ALUMINUM 


K = 18  ksi  in.^ 
max 

K . =9  ksi  in. 

min 


Equilibrium  Superdislocation  Position  ( 10~ 3 inches) 


Load  Cycle 


ho 


Superdislocation  Strengths  and  Crack-Growth  Increment 


Load  Cycle 


7 


Linearized  versions  of  the  dislocation  interaction  terms  can  be  deduced 
by  assuming  that  both  dislocations  are  on  the  same  slip  plane  and  that  they 
occupy  positions  such  that  z » 23  » lz  . Hence, 


g32  + k3  2 = - (1-  cos  4 0)  - 7 r (1  + cos  0)  sin°  0 

il  2 Ux  *2)2 

and 

3 1 2 * - l 

Ssi  + ^21  = o'  “3  (1  +cos  9)  sins  0+2  ~T  (1  - cos  4 9)  . 

1 


For  convenience  in  the  following,  the  specific  values  for  all  the  linearized 
functions  for  0 = cos-1  (1/3)  can  be  written  out.  These  are 


fx  = 6.532  (a/j^ 

§n  + ^li  = 

h31  = 0.2722  (a/l-i  + 0.314 

§12  + bl  2 = " ~ 


f2  = 6.532  (o/j J3)% 
§22  ^22  = " 1/^2 


h2  = 0.2722  (a/2s)2 
1 


(-Ms)1 


128  JL_ 

81  X, 


+ 0.314 


§21  + k21 


16  128  fs 

9 23s  + 81  V ‘ 


A comparison  showing  the  accuracy  of  the  latter  two  functions  is  shown  in 
Table  4.3. 


The  basic  premise  of  the  steady-state  growth  model  is  that  the  strength 
and  position  of  the  super-superdislocation  are  completely  determined  by  the 
conditions  of  force  equilibrium  at  maximum  and  minimum  load  in  the  cycle.  It 
is  assumed  that  the  resistance  to  outward  motion  just  balances  the  stress  act- 
ing to  push  the  superdislocation  out  at  maximum  load.  At  the  minimum  load,  the 
resistance  just  balances  the  inward  driving  stress.  The  appropriate  relations 
therefore  are 


T , 
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o h + 
max 


G 

TT  (<  + 1) 


b {g  + k] 


-T 
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a . 

min 


h 


+ 


G 

TT  (*+  1) 


b { g + k]  . 


(4.26) 
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TABLE  4-3.  COMPARISON  OF  EXACT  AND  APPROXIMATE  FUNCTIONS 


Exact 


Approximate 


o-3 

o"3 

-a 

0“3 

-a 

-l 

-l 

-l 


10“* 

10“5 

10“4 

10“3 
10“5 
10  “4 

10“3 
10“5 
10  “4 


-57799 
-1937A 
- 7386 


-17950 

- 5784 

- 1939 


57.79 

193.71 

738.48 


1.79 

5.77 

19.36 


-57798 

57.80 

-19358 

193.58 

- 7202 

720.21 

17936 

1.79 

5780 

5.78 

1936 

19.36 

- 5684 


- 1808 


- 5638 


Using  Che  linearized  version  of  the  functions  appearing  in  these  equations  and 
taking  the  case  9 = cos"1 (1/3),  then  gives 


(4.27) 


Equations  (4.27)  can  be  readily  solved  to  obtain  b and  l.  The  results  are 

. Vl  -s 


l = 


T . 
1 


54  a a . 9 

max  - min 


(4.28) 


where  l is  assumed  to  be  positive  and 


, _ 7r(<+m  j [7  2 ^ 8*1  } 

b G )CTmaxLl27y  U)  + 9 J " Ti  j ■ 


(4.29) 


Equivalently,  using  yi g and  /*a  and  taking  x.  = Y//3, 


max  max 
Equation  (4.28)  can  be  written 


, /K  - K , 

1 J max  min 

3| 

1 f-L? 

O (J  . \) 

max  - mm  1 

18tt  \ Y 

1 

1 V 27/ 

1 [( 

-2 


(4.30) 


whereupon  Equation  (4.29)  becomes 


b = ~ 


§ & 4K{K  -f 
27  YG  ) max  2 


1 -(A' 

'27' 

.t  0 

2 max 
Y 

1 

1 

1-  (-L)  Afl. 

V277  Y 


-1 


, (4.31) 


where  Act  = a - ct  . , and  AK  = K - K . 


Now,  assuming  Che  crack-growch  rate  to  be  as  discussed  above,  then  it  follows 


that 


— = r b 
dN  27/3  ^ 


cos 


e . 


(4.32) 


where  C is  a (dimensionless)  number  to  be  determined  empirically  as  follows. 

Consider  plane-strain  conditions  where  k + 1 = 4(l-v)  and,  for 
simplicity,  assume  that  <^m3X  « Y.  Then,  substituting  Equation  (4.31)  into 
Equation  (4.32)  gives 


da  _ 4 _ 

dN  “ 27/3  * 


n - v 

V YG 


) lKmax 


(4.33) 


Note  that  the  material's  properties  enter  this  relation  only  through  the  param- 
eter (l-v)/YG.  Typical  values  are  given  in  Table  4.4. 


TABLE  4.4.  VALUES  OF  THE  MECHANICAL  PROPERTIES  USED  IN  THE 


COMPUTATIONS 

OF  FATIGUE 

-CRACK- 

GROWTH  RATES 

1 - V 

Material 

G (ksi) 

V 

Y (ksi) 

GY 

Ti-6A1-4V 

6000 

0.33 

120 

0.931  x 10" ^ 

2024-T3A1 

3750 

0.34 

52.5 

3.35  x 10"s 

The  empirical  parameter  £ can  be  estimated  by  the  results  obtained  using 
the  complete  cycle-by-cycle  computational  model.  The  basis  for  this  estimate 
is  given  in  the  following  section  of  this  report.  The  result  is 

£ = 0.200  (dimensionless)  . (4.34) 

Hence,  the  approximate  version  of  the  steady-state  crack-growth  relation  that 
is  the  objective  of  this  section  can  be  written 


d; 

dN 


,018  ( 


n - v\  r 


V YG 


) ;K 


AK  - \ 
max  2 


(AK)2 


(4.35) 


where,  for  definiteness,  the  parameters  appearing  in  this  relation  can  be 
defined  explicitly  as  follows: 
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G 

Y 


= shear  modulus  (ksi  or  MN/m2  ) 

= tensile  yield  stress  (ksi  or  MN/nf ) 


v 


K 

max 

K . 
mm 

AK 


Poisson's  ratio  (dimensionless) 

stress-intensity  factor  at  maximum  load  (ksi  in. 2 or  MNm  2) 

\ «.-2. 

stress-intensity  factor  at  minimum  load  (ksi  in. or  MNm  c) 

K - K . 
max  min 


Using  numerical  values  in  the  English  system  of  units  in  Equation  (4.35)  will 
give  crack-growth  rates  in  units  of  in. /cycle,  the  metric  units  noted  will 
give  growth  rates  in  M/ cycle. 

An  alternative  way  of  writing  Equation  (4.35)  is  by  introducing  the  stress 
ratio  parameter,  R.  This  is 


(4.36) 


This  parameter  can  be  used  to  replace  K since  from  Equation  (4.36),  K 

^ v max  M max 

AK/(1-R).  Equation  (4.35)  can  then  be  written  as  a function  of  R and  AK  as 


1^'J  (I±^,  (AK)2 
, YG  / \1  - R/  U ' ' 


(4.37) 


This  alternative  form  of  the  steady-state  fatigue-crack-propagation-rate  rela- 


tion arising  from  the  super-superdislocation  approach  is  more  useful  for  the 
cases  where  R is  a constant. 


Discussion  of  the  Steady-State  Crack 
Growth-Rate  Model  Obtained  With  the 
Super-Superdis location  Approach 


It  may  be  useful  to  consider  briefly  some  qualitative  features  of  Equation 
(4.35).  The  first  point  is  in  regard  to  the  form  of  the  equation.  It  can  be 
seen  that  the  AK  dependence  corresponds  to  a "power-law"  growth  rate  with  the 
power  2.  It  is  worth  remarking  that,  at  this  stage  of  the  development  of  the 
model,  this  result  could  be  expected  with  absolute  certainty  purely  from 
dimensional  analysis.  That  is,  there  are  only  two  length  parameters  that  have 
been  allowed  in  the  linearized  version  of  the  problem:  Kf  and  b.  Hence,  the 

resulting  expression  for  da/dN  must  involve  K2  only. 


Most  empirically  derived  steady-state  fatigue-crack-growth  relations 
incorporate  an  exponent  in  a power-law  type  of  relation  that  is  more  nearly 
equal  to  3 than  to  2.  For  example,  the  extensive  results  on  2024-T3  aluminum 
reported  by  von  Euw,  et  al  [7],  are  very  well  represented  by  the  simple  rela- 
tion (English  units) 


da 

dN 


4.32  x lO"9 


OK)3 
1 - R 


(4.38) 


This  indicates  that  at  least  one  additional  parameter  with  a length  dimension 
must  somehow  enter  the  problem.  Anticipating  that  this  parameter  could  be 
connected  with  the  microstructure  of  the  material,  a survey  of  the  literature 
on  the  micromechanical  aspects  of  fatigue-crack  growth  was  performed  as  part 
of  this  project.  The  preliminary  results  of  this  study  are  given  in  Appendix 
C of  this  report. 

The  second  point  is  connected  with  the  empirical  parameter  £ introduced 
in  Equation  (4.32)  to  account  for  the  fact  that  only  a small  part  of  the 
residual  plasticity  represented  by  the  super-superdislocation  is  associated 
within  the  crack  tip  itself.  That  is,  £ represents  the  "sliding  off"  of  the 
dislocations  on  slip  planes  emanating  from  the  crack  tip  that  can  be  identified 
with  the  forward  thrust  of  the  crack  during  a given  load  cycle.  As  discussed 
above,  this  length  is  exactly  equal  to  the  strength  of  the  superdislocation. 
Hence,  the  forward  component  of  crack  growth  is  simply  b cos  e.  In  the  super- 
superdislocation model,  the  parameter  £ can  be  given  the  physical  interpreta- 
tion of  being  the  ratio  of  the  strength  of  the  superdislocation  on  the  crack 
tip  slip  plane  to  the  total  strength  of  the  super-superdislocation. 

A rough  estimate  of  the  magnitude  of  £ can  be  made  by  examining  the 

results  of  the  cycle-by-cycle  computational  model.  Assuming  for  simplicity 

that  the  strength  of  the  nascent  superdislocation  is  the  same  as  for  those 

created  in  each  previous  load  cycle,  the  number  of  load  cycles  that  contribute 

to  the  steady-state  process  is  £-1  . Consequently,  £_1  load  cycles  should  be 

necessary  to  establish  the  steady  state.  Now,  as  exemplified  by  the  results 

shown  in  Figure  5.3  (see  following  section),  a steady-state  growth  rate  was 

in  fact  achieved  in  about  five  cycles  of  loading.  This  was  true  for  a range 

of  K and  K . values.  Hence,  £ can  be  estimated  to  be  0.2. 

max  min 

To  be  sure,  this  comparison  is  very  approximate.  Nevertheless,  it  does 
provide  some  reassurance  that  the  good  agreement  with  experiment  described  in 


k 


Che  following  is  not  contrived  but  has,  instead,  a fundamentally  sould  basis. 
It  also  suggests  that  a direct  computation  of  £ can  be  made  from  the  complete 
cycle-to-cycle  model. 

The  third  point  is  in  the  appearance  of  purely  well-known  mechanical 
properties  of  the  material  in  Equation*  (4. 35 ) . Aside  from  the  possible  slight 
dependence  of  £ on  the  material,  this  means  that  the  growth  rate  "law"  derived 
from  the  super-superdislocation  approach  can  be  easily  applied  to  any  material 
without  the  necessity  of  performing  special  tests  to  determine  material- 
dependent  empirical  parameters  required  in  the  relation.  While  exhaustive 
comparisons  with  experiments  have  not  been  performed  in  this  work,  some 
verification  is  offered  in  the  next  section  where  it  is  shown  that  the  model 
works  equally  well  for  both  2024-T3  aluminum  and  Ti-6A1-4V  titanium  alloy. 

Finally,  while  not  in  precise  agreement  with  all  fatigue-crack-growth 
data.  Equation  (4.35)  does  reflect  a number  of  experimental  observations.  In 
particular,  for  the  most  common  experimental  conditions,  AK  = constant.  Equa- 
tion (4.35)  would  predict  with  reasonable  accuracy  a linear  dependence  of  the 
growth  rate  on  Of  most  importance  in  regard  to  the  eventual  use  of  the 

superdislocation  concept  in  a complete  fatigue-crack-growth  model  is  that  no 

arbitrary  disposable  parameters  need  to  be  introduced  into  the  model  to  obtain 
d2 

a prediction  of  — . As  exemplified  by  Equation  (4.35),  the  only  parameters 
that  enter  into  the  analysis  are  the  basic  mechanical  properties  G,  Y,  and  v. 
This  feature,  together  with  the  ability  of  the  model  to  admit  arbitrary  varia- 
tions in  the  load  level  from  cycle  to  cycle,  make  this  approach  unique. 

5.  COMPARISONS  BETWEEN  ANALYSIS  AND  EXPERIMENT 
Comparison  With  Alternative  Analyses 


Two  features  that  likely  control  the  rate  of  fatigue-crack  growth  are 
(1)  the  size  of  the  intensely  deformed  plastic  region  surrounding  the  crack 
tip  that  is  created  in  each  load  cycle  and  (2)  the  extent  to  which  the 
plastically  deformed  regions  created  in  different  load  cycles  interact  with 
each  other.  The  first  of  these  can  be  associated  with  the  CTOD  calculated 
with  a strip-yield  model.  It  might  be  noted  that  strip-yield  models  of  the 
Dugdale  type  also  offer  predictions  of  the  COD.  However,  in  the  Dugdale  model, 
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plastic  zone  interactions  do  not  take  place  in  a realistic  way.  At  the  other 
extreme,  a finite  element  or  other  two-dimensional  elastic-plastic  computational 
scheme  applied  to  fatigue-crack  growth  requires  so  much  computer  time  as  to  be 
completely  impractical  for  more  than  a few  load  cycles.  Thus,  the  analytical 
approach  taken  in  this  work  is  the  simplest  that  takes  direct  account  of  the 
essential  features  of  fatigue-crack  growth  and,  at  the  same  time,  permits  a 
large  number  of  arbitrary  load  cycles  to  be  studied. 

The  analytical  approach  taken  in  this  work  rests  on  two  fundamental  assump- 
tions. The  first  is  that  the  CTOD  arising  from  the  intense  plastic  deformation 
attending  the  crack  tip  is  the  parameter  that  governs  crack  growth.  The  second 
is  that  the  essential  interaction  between  the  plastic  deformation  created  in 
different  load  cycles  can  be  represented  by  an  inclined  strip-yield  superdislo- 
cation model.  The  approach  taken  here  uses  what  is  perhaps  the  simplest 
mathematical  model  in  which  those  considerations  (i.e.,  CTOD  and  distinguish- 
able plastic  zones  in  each  load  cycle)  can  be  taken  into  account  in  a rigorous 
fashion. 

In  the  superdislocation  approach  as  applied  to  a fatigue  situation,  the 
plastically  deformed  regions  surrounding  the  crack  tip  are  represented  by 
single  pairs  of  dislocations.  These  are  considered  to  "slide  off"  from  the 
crack  tip  on  slip  planes  inclined  at  a fixed  angle,  0,  from  the  crack  plane. 

This  is  shown  in  Figure  5.1.  It  can  be  seen  that  the  plastic  deformation 
resulting  from  a series  of  cyclic  load  applications  is  discrete,  giving  a 
"jogged"  appearance  to  the  crack  surfaces.  The  "sliding  off"  at  the  crack  tip 
gives  rise  to  a CTOD,  5,  and  suggests  a value  for  the  increment  of  crack 
advance  per  load  cycle,  Afl . It  is  worth  reemphasizing  that  crack  closure,  if 
it  is  to  occur,  will  occur  naturally  in  this  model  and  need  not  be  imposed 
arbitrarily.  Conceptually,  the  analytical  approach  taken  here  stands  somewhere 
between  the  approximate  semiempirical  treatments  of  Wheeler  [8]  and  of  Willen- 
borg,  et  al  [9],  and  the  more  precise,  but  cumbersome,  finite-element  elastic- 
plastic  computation  being  developed  by  Newman  [10].  To  put  the  model  into 
better  perspective,  a brief  description  of  these  competing  approaches  to 
fatigue-crack  propagation  under  variable-amplitude  cyclic  loading  may  be  in 
order.  As  background,  the  four  principal  models  of  crack-tip  plasticity  are 
shown  in  Figure  5.2.  These  are  listed  in  order  of  both  increasing  complexity 
and  accuracy.  (It  may  be  worth  noticing  that  only  the  strip-yield  models,  (b) 
and  (c),  give  an  unambiguous  prediction  of  the  CTOD.)  Four  classes  of 
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(a)  Modified  Linear  Elastic-Fracture  Mechanics 


Linear  elasticity  solutions  for  partly  loaded 
crack  with  surface  tractions  taken  equal  to 
tensile  yield  stress  — singularity  is  cancelled. 


(b)  Dugdale  (BCS)  Strip-Yield  Model 


Dislocation  model  with  slip  confined  to  pair  of 
slip  planes  emanating  from  the  crack  tip.  Singu 
larity  is  cancelled.  Superdislocation  model  is 
simple  variation  of  this. 


(c)  Bilby-Swinden  Strip-Yield  Model 


Typical  "butterfly"  yield  zone  at  tip  of  crack 
under  plane-strain  conditions  for  elastic/ 
perfectly  plastic  material.  Singularity  not 
cancelled  in  linear  theory.  Large  strain  theory 
shows  existence  of  "intensely  strained  region" 
at  crack  tip  where  strains  — unity. 


(d)  Plane  Strain  Elastic-Plastic  Model 


FIGURE  5.2.  COMPARISON  OF  CRACK-TIP  PLASTICITY 
MODELS  USED  IN  FATIGUE  AND  FRACTURE 
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techniques  for  fatigue-crack-growth  prediction  can  be  identified  on  the  basis 
of  the  categorization  given  in  Figure  5.2.  These  are  listed  and  described  in 
Table  5.1.  Note  that  these  can  similarly  be  numbered  in  order  of  increasing 
accuracy  and  complexity. 

Other  models  which  might  be  mentioned  are  those  of  Yokobori  and  Yoshida  [15] 
and  of  Neumann  [16].  While  these  have  some  resemblance  to  the  present  work, 
they  have  not  exploited  the  potential  of  this  approach  to  nearly  the  extent 
given  in  this  report. 

Comparison  With  Experimental  Results 


f ! 


\ 
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Qualitative  results  that  can  be  compared  with  experimental  results  and 
with  analytical  results  obtained  from  more  precise  approaches  (e.g.,  finite- 
element  methods)  have  been  obtained  in  two  areas 

(1)  Equilibrium  under  a single  load  cycle 

(2)  Crack  growth  under  constant  AK  cycling. 

Using  the  notation  developed  in  the  preceding  section,  the  pertinent  results 
for  "small-scale  yielding"  are  found  to  be 


i = .01354  -rr  , 
Y* 


and 


a = 0.4763 


1 - vs  K2 


(5.1) 


(5.2) 


where  9 has  been  selected  as  the  angle  of  maximum  shear  stress  to  maximize 
the  extent  of  yielding  (i.e.,  9 = cos-1  (1/3)  = 70.3°)  and  the  internal  fric- 
tion stress  has  been  taken  in  Y//7.  Now,  while  these  results  are  not  important 
in  themselves,  they  provide  a check  on  the  fidelity  of  the  model.  In  partic- 
ular, many  different  analytical  approaches  have  been  employed  to  investigate 
the  small-scale  yielding  regime  and  all  are  in  qualitative  agreement  with 
Equation  (5.2).  Quantitatively,  as  the  comparison  given  in  Table  4.1  clearly 
shows,  the  models  also  give  quite  acceptable  results.  Again,  these  findings 
do  not  have  any  direct  bearing  on  the  application  of  the  model  to  fatigue-crack 
growth  except  as  a necessary  basic  check. 

The  extension  of  the  model  to  multicycle  loading  has  been  described  in  the 
preceding  section  of  this  report.  It  need  only  be  pointed  out  here  that  to 
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perform  a fatigue-crack-growth  propagation,  all  that  is  needed  is  a specific 
cyclic  loading  sequence  and  three  material  parameters:  the  shear  modulus,  G; 

Poisson's  ratio,  v;  and  the  tensile  yield  stress,  Y.  Computations  have  been 
performed  for  both  an  aluminum  and  a titanium  alloy  using  the  property  values 
given  in  Table  5.2.  It  is  believed  that  the  differences  between  these  two 
sets  of  values  are  large  enough  to  test  the  ability  of  the  model  to  distinguish 
between  the  two  materials. 

TABLE  5.2.  MATERIAL  PROPERTIES  USED  IN  COMPUTATION 
OF  FATIGUE-CRACK  PROPAGATION 


Material 

G (ksi) 

V 

Y (ksi) 

2024-T3  Aluminum 

3750 

0.34 

52.5 

Ti-6A1-4V  Titanium 

6000 

0.33 

120 

Example  results  are  shown  in  Figure  5.3  for  the  titanium  alloy.  Plotted 

in  Figure  5.3  are  the  crack-growth  increments  for  each  load  application  for 

three  different  values  of  IK  = K - K . for  the  special  case  where  R = 

max  mm  r 

K . /K  = 0.5.  The  most  important  features  of  the  results  shown  in  Figure 
min  max  r 

5.3  are,  first,  that  an  "accelerated"  growth  is  evident  in  the  first  few  load 
cycles  and,  second,  that  the  growth  rate  levels  off  to  a constant  rate  charac- 
teristic of  the  AK  value  applied.  These  two  features  are  in  good  qualitative 
agreement  with  observations  of  fatigue-crack  growth.  A quantitative  comparison 
based  on  the  steady-state  growth  rates  for  both  the  titanium  alloy  and  the 
aluminum  alloy  are  shown  in  Figures  5.4  and  5.5,  respectively.  On  the  basis 
of  the  agreement  evident  in  these  two  figures,  it  is  believed  that  the  model 
has  proven  to  be  more  than  adequate,  particularly  in  view  of  the  fact  that  no 
arbitrary  factors  have  been  introduced  to  adjust  these  results. 

With  regard  to  the  crack-closure  phenomenon,  attention  can  be  given  to 
the  calculated  crack-face  displacements.  Example  results  are  shown  for 
titanium  in  Figure  5.6  and  for  aluminum  in  Figure  5.7.  Here,  the  crack-open- 
ing displacements  are  shown  at  maximum  and  minimum  load  in  the  twelfth  load 
cycle.  The  interval  shown  is  that  through  which  the  crack  has  grown  since  the 
initiation  of  growth  at  the  first  load  cycle.  It  can  be  seen  that  the  crack 
faces  are  in  fact  closing  upon  each  other  during  the  minimum  load. 
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Curve  C 


Curve  B 


Curve  A 


Number  of  Load  Cycles 


COMPUTATIONAL  RESULTS  FOR  FATIGUE-CRACK  GROWTH  IN  TITANIUM 
USING  INCLINED  STRIP-DISLOCATION  MODEL  FOR  R = 0.5 


FIGURE  5.3 


Curve 
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Model  Parameters 

G = 6000  ksi 
v = 033 
Y = 120  ksi 


Scatter  band  of 
experimental  results 


Model  Parameters 

G = 3750  ksi 
u - 0.34 
Y = 52.5  ksi 


Computational  result 


FIGURE  5.5 


COMPARISON  OF  COMPUTATIONAL  AND  EXPERIMENTAL  RESULTS  IN 
2024-T3  ALUMINUM  FOR  STEADY-STATE  CONDITIONS  WITH  R = 0 
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On  the  negative  side,  there  are  certain  predictions  of  the  present  model 

that  are  not  in  such  good  agreement  with  experiment.  For  example,  while  the 

model  predicts  an  "R-effect",  it  is  contradicted  by  the  data.  That  is,  while 

the  growth  rates  for  uniform  cyclic  loading  with  AK  = constant  tend  to  increase 

with  R = K . /K  , the  model  predicts  a decrease.  This  is  shown  in  Figure 
mm  max 

5.8  for  aluminum. 

The  trend  for  the  K . =0  calculations  to  overestimate  the  growth  rate, 

mm 

seen  in  Figure  5.8,  is  also  shown  in  Figure  5.9.  In  the  latter  figure,  the 

measured  crack-growth  rates  for  aluminum  as  a function  of  K for  K . =0  can 

max  mm 

be  seen  to  be  considerably  under  the  predicted  values.  However,  the  correct 
general  trend  does  seem  to  be  given  by  the  model.  Subsequent  work  is  to 
improve  upon  the  deficiencies  in  the  predictions  of  the  current  model. 

The  Effect  of  Mechanical  Properties 


As  has  been  emphasized  throughout  this  report,  the  computational  model 
depends  only  on  the  ordinary  mechanical  properties  of  the  material  being 
simulated:  G,  Y,  and  v.  It  is  of  interest  to  determine  the  relative 

importance  of  each  of  these  on  the  fatigue-crack-growth  rates.  From  consider- 
ation of  the  basic  equations  of  the  model,  Poisson's  ratio  can  be  seen  to  have 
little  effect.  So,  only  G and  Y need  to  be  varied.  This  has  been  done  by 
taking  the  particular  values  used  for  the  aluminum  and  titanium  simulations 
and  interchanging  them.  The  results  are  as  shown  in  Table  5.3. 


TABLE  5.3.  COMPARISON  OF  STEADY-STATE  CRACK-GROWTH  RATES  SHOWING 

THE  EFFECT  OF  THE  MECHANICAL  PROPERTIES  OF  THE  MATERIAL 


K = 20,  K . =10 

max  mm 


v = 0.335 

/d2\ 

VdijJ  un^ts  * in. /cycle 


Y = 52.5  ksi 


Y = 120.0  ksi 


•vJ 


G = 3750  ksi 
G = 6000  ksi 


9.61 

5.91 


4.10 

2.48 


Steady-state  crack  growth  rate 
predicted  by  cycle-by-cycle 
analysis  model 


Denotes  experimental 
data  of  Von  Euw,  et  al, 
ASTM  STP  513  (1972) 
pp  230-259. 


Kmax>  ksHn- 

FIGURE  5.8.  COMPARISON  OF  PREDICTED  AND  OBSERVED  FATIGUE-CRACK-PROPAGATION 
RATES  IN  2024-T3  ALUMINUM  FOR  AK  = 9 KSI-IN.^ 
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Steady -state  crack  growth  rate 
predicted  by  cycle-by-cycle  ^ 
analysis  model 


O Denotes  experimental 
data  of  Von  Euw,  et  at, 
ASTM  STP  513  (1972) 
p 230. 


COMPARISON  OF  PREDICTED  AND  OBSERVED  FATIGUE-CRACK- PROPAGATION 
RATES  IN  2024-T3  ALUMINUM  FOR  K^  = 0 


FIGURE  5.9 


The  results  in  Table  5.3  demonstrate  a (GY)-1  dependence  of  the  fatigue- 
crack-growth  rates  in  the  cycle-by-cvcle  computations.  This  was,  of  course, 
already  seen  in  the  closed-form  expressions  found  using  the  super-superdislo- 
cation model  in  the  preceding  section.  However,  because  further  modifications 
and  extensions  of  the  model  are  planned  to  correct  the  R-effect  discrepancy 
in  the  present  model,  it  is  conceivable  that  this  dependence  might  also  be 
altered  subsequently. 


A Heuristic  Crack-Growth  Retardation  Calcula 
tion  Using  the  Super-Superdislocation  Model 


The  super-superdislocation  concept  introduced  in  Section  4 can  be  used  to 
demonstrate  how  the  model  can  be  made  to  reflect  crack-growth  retardation 
following  an  overload.  A striking  example  of  this  phenomenon  obtained  by 
von  Euw,  et  al  [7],  in  experiments  on  2024-T3  aluminum  is  shown  in  Figure  5.10. 
The  super-superdislocation  model  can  be  applied  in  a heuristic  fashion  by 
making  two  simplifying  assumptions.  First,  that  the  residual  plastic  deforma- 
tion contributed  by  all  the  load  cycles  following  the  overload  can  be  neglected 
in  comparison  to  the  overload  plasticity.  Hence,  the  latter  can  be  represented 
by  a single  fixed  super-superdislocation.  The  second  assumption  has  to  do  with 
the  effect  of  the  super-superdislocation  on  a growing  crack  tip.  Specifically, 
it  is  assumed  that  the  maximum  effect  of  the  super-superdislocation  is  felt 
on  the  crack  plane  at  the  point  directly  under  it.  Its  effect  is  assumed  to 
diminish  thereafter  and  become  zero  at  the  point  where  a normal  to  the  slip 
plane  drawn  from  the  super-superdislocation  position  intersects  the  crack 
plane.  These  two  points  are  indicated  in  Figure  5.11. 

As  can  be  seen  in  Figure  5.11,  the  crack  length  at  which  the  maximum 
effect  occurs--and  where  the  minimum  growth  rate  is  expected--is  at  a distance 
H cos  0 from  the  crack  length  at  the  imposition  of  the  overload.  Similarly, 
the  original  growth  rate  should  be  reached  after  a growth  increase  of  l/cos  0. 
Now,  taking  the  linearized  equations  for  the  single  dislocation  pair  and 
solving  for  the  position  gives 


AK  = I8  ksi  in.l/2,  100% 
peak  overload 


-0.04  0 0.04  0.08  0.12  0.16  0.20 

Distance  From  the  Peak  Load,  a,  in. 

FIGURE  5.10.  RETARDED  FATIGUE- CRACK- GROWTH  RATE  IN 
2024-T3  ALUMINUM  ALLOY  FOLLOWING  THE 
APPLICATION  OF  A SINGLE  OVERLOAD^] 


Super-superdislocation  slip  plane 

Super-superdislocation  representing  an 
overload  cycle 


Notice  that  if  = Y sin  9 cos  0,  this  has  the  very  attractive  feature  of 
giving  an  infinite  extent  of  yielding  and  infinite  extent  of  yielding  and 
infinite  COD  at  a = Y,  just  as  in  the  Dugdale  model.  Accepting  this  for  the 
moment  as  being  more  physically  plausible  than  the  choice  used  previously, 
the  potential  of  the  model  to  reproduce  some  features  of  growth-rate  retarda- 
tion can  be  examined. 

Taking  a - K ) , t.  = Y sin  9 cos  9 and  inserting  0 = cos-1  (1/3), 

OL  O X 

the  above  expression  reduces  to 


£ * -0*57  (- ff  (l -1 


For  K = 36  ksi  in. 2 (the  overload)  and  Y = 52.5  ksi  (aluminum)  and  taking  20  = 
1 inch,  this  gives  l = .057  inch.  Now,  the  increase  in  crack  length  such  that 
the  crack  tip  is  directly  under  the  dislocation  produced  in  the  overload  cycle 
is  = l cos  9 = £/3  = .019  inch.  In  comparing  with  Figure  5.10,  it  can  be 
seen  to  be  about  equal  to  the  distance  from  the  peak  load  at  which  the  minimum 
growth  rate  occurs.  The  increase  in  crack  length  such  that  the  center  line  of 
the  dislocation  intersects  the  crack  tip  is  = U cos  9 =3 1,  or  0.172  inch. 
Again  referring  to  Figure  5.10,  this  is  just  about  the  distance  from  the  peak 
load  at  which  the  influence  of  the  peak  load  vanishes  and  the  steady-state 
growth  rate  is  recovered. 

Comparisons  can  also  be  made  with  the  data  discussed  in  Section  3 of  this 
report  for  the  effect  of  peak  overloads  in  Ti-6A1-4V  titanium  alloy.  For  this 
purpose,  Equation  (5.3)  can  be  used  to  express  the  length  of  crack  affected  by 
the  overload  ir.  the  form 


CL  = 0.1371 


f Y 


- -a 
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oL 
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(5.4) 


which  reveals  that  the  dependence  is  not  simply  a function  of  the  stress- 
intensity  factor. 

Figure  5.12  shows  a comparison  between  the  experimental  data  and  the 
prediction  of  Equation  (5.4)  for  various  assumed  values  of  c0 . It  must  be 
pointed  out  that  the  agreement  is  deceptive  because  the  actual  data  were 
generated  with  generally  larger  values  of  £0 . Admittedly,  there  is  consider- 
able scatter  evident;  however,  the  comparison  being  sought  is  a qualitative 
correlation.  The  actual  30  values  used  were  on  the  order  of  1 inch  and  larger 
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FIGURE  5.12 


DISTANCE  OF  CRACK  PROPAGATION  FOLLOWING  AN  OVERLOAD  BEFORE 
STEADY -STATE  GROWTH  RATE  IS  RECOVERED  IN  Ti-6A1-4V 


for  experimental  sensitivity.  The  lower  values  shown  on  the  computed  lines 
are  used  to  display  the  predicted  trend  in  the  more  design-sensitive  range. 
Nevertheless,  in  view  of  the  simplicity  of  the  present  calculation,  the 
qualitative  agreement  is  quite  instructive.  It  can  be  concluded  that  the 
model  does  have  the  potential  to  accurately  predict  crack-growth  retardation. 


Closure 
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In  this  section  of  the  report,  computations  have  been  performed  for  two 
different  materials:  2024-T3  aluminum  and  Ti-6A1-4V  titanium.  For  loading 

corresponding  to  AK  being  constant,  the  results  show  that  a relatively  constant 
fatigue-crack-growth  rate  is  quickly  achieved;  e.g.,  within  about  ten  load 
cycles.  Although  some  discrepancies  exist,  the  comparisons  between  these 
steady-state  growth  rates  and  experimental  results  for  these  materials  show 
quite  good  agreement. 

For  reasons  that  are  not  entirely  clear  at  this  point,  computations  made 
for  the  special  case  where  R = K . /K  =0.5  are  virtually  right  in  the 
middle  of  the  scatter  band  of  the  experimental  data.  This  is  true  for  both 
the  aluminum  and  the  titanium  materials  for  a range  of  AK  values.  On  the 
negative  side,  while  the  calculations  do  show  both  an  "R  effect”  and  a load 
interaction  effect  in  nonconstant  amplitude  cycling,  in  neither  case  are  the 
results  satisfactory  as  yet.  These  are  probably  connected  so  that  the  remedy 
for  one  will  likely  make  a substantial  improvement  in  the  other,  however.  A 
number  of  options  are  open  in  which  improvement  of  the  model  can  be  made. 

These  are  discussed  in  the  next  section. 

6.  CONCLUSIONS  AND  RECOMMENDATIONS 
FOR  FURTHER  RESEARCH 


The  specific  objective  of  this  work  was  to  apply  a simple  model  of  plastic 
relaxation  at  a crack  tip  to  the  fatigue  problem.  The  model  allows  an  estimate 
of  the  plastic  flow  at  each  stage  in  a reasonably  simple  way  using  the  "super- 
dislocation"  idea  which  further  allows  a look  at  cycle-bv-cycle  growth  with 
residual  plasticity  left  from  previous  cycles  still  playing  a role.  Thus, 
a sequential  calculation  is  made  in  a logically  self-consistent  manner  with 
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due  attention  being  paid  to  the  equilibrium  of  each  of  the  superdislocations. 
Perhaps  the  main  criticism  that  could  be  leveled  at  this  approach  is  that  the 
superdislocation  is  an  oversimplification  of  the  plastic  field  which  may  limit 
the  chances  of  the  present  model  exhibiting  the  trends  usually  observed  in 
experiments  on  fatigue.  In  fact,  as  pointed  out  earlier  in  this  report,  the 
present  results  do  not  in  all  respects  duplicate  those  to  be  obtained  by 
experiment.  However,  this  is  due  to  the  details  contained  in  the  present 
model  rather  than  to  the  overall  philosophy.  In  this  section,  some  of  the  main 
consequences  of  the  results  are  discussed  and  recommendations  for  further  work 
to  overcome  the  objectionable  features  are  given. 

In  the  first  place,  with  the  complete  model,  a strictly  logical  evalua- 
tion is  conducted  with  the  results  being  a direct  consequence  of  the  growth 
criterion  and  the  variety  of  stress  interactions  occurring  between  the  dislo- 
cations and  the  crack.  As  one  highly  significant  result,  the  cycle-by-cycle 
calculations  do  predict  a steady-state  growth  rate.  Before  discussing  possible 
changes  that  could  be  made  in  the  growth  conditions,  consider  the  computational 
results  in  terms  of  a simplified  view  of  the  steady-state  situation  as  applied 
to  the  complete  model.  The  way  that  the  growth  is  specified  on  the  crack  at 
the  a . load,  the  crack  closes  at  the  crack  tip  after  advance  (4  = 0)  and 
the  dislocation  on  the  current  slip  band  changes  its  strength  and/or  position 
to  accommodate  the  condition  as  well  as  being  in  equilibrium.  The  net  effect 
is  that  when  ^max  is  applied  again,  it  is  acting  on  a residual  stress  due  to 

the  dislocations  left  behind.  These  residual  stresses  iust  balanced  a . . 

J min 

Hence,  on  the  current  crack,  an  effective  stress  of  (a  - a . ) is  acting. 

’ max  min 

It  is  then  expected  from  the  linearized  results  that  6 = 0.44  K2 / EY . This 

result  for  crack  advance  is  proportional  to  iK  = K - K . = /tto  (a 

r r max  mm  max 

<7  . ) . This  neglects,  of  course,  any  effect  that  neighboring  dislocations 
min 

have  on  the  current  one  being  emitted  from  the  crack  tip;  therefore,  the  result 
can  only  be  approximate.  The  actual  computational  results  are  quite  close  to 
the  results  suggested  by  the  above  argument.  Hence,  the  numerical  results 
for  the  complete  model  demonstrate  that  the  progressive  development  and 
dispersion  dislocations  will  in  fact  lead  to  a steady-state  growth  rate  which 
is  approximately  proportional  to  (AK)2  even  when  all  the  dislocation  stress 
fields  are  correctly  taken  into  account. 

Previous  models  appearing  in  the  literature  have  never  shown  this  with  a 
self  consistent  logical  scheme  in  which  the  process  of  producing  dislocations, 
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etc.,  is  followed  through  in  a sequential  way.  It  might  be  noted  that  the 
model  does  demonstrate  a deviation  from  the  (AK)2  law  which  depends  precisely 
on  the  dislocation  interactions  (our  model  of  the  plastic  field).  However, 
the  precise  form  of  the  complete  numerical  results  do  not  show  the  desired 
experimental  trend.  So,  some  modifications  are  obviously  in  order. 

As  discussed  in  the  preceding  section  of  this  report,  there  are  many 
aspects  of  the  present  model  which  can  certainly  be  improved  upon.  By 
exploring  the  options  that  are  open  to  improve  the  present  shortcomings  and 
incorporating  the  satisfactory  candidates  to  build  upon  the  presently  success- 
ful features,  it  is  believed  that  these  can  be  overcome.  Consequently,  the 
model  can  be  used  as  the  basis  of  the  computational  approach  that  is  the 
objective  of  this  work.  Some  of  the  possibilities  are  as  follows: 

(1)  The  singularity-canceling  equation  is  currently  imposed  at 
the  maximum  as  well  as  at  the  minimum  load  level  in  each 
load  cycle.  In  addition  to  being  physically  unreasonable, 
this  unduly  constrains  the  crack-tip  deformation  and,  hence, 
leads  to  incorrect  growth-rate  predictions.  As  one  example,  a 
threshold  AK  level  is  precluded.  What  is  essentially  needed 
is  a mathematical  relation  that  can  be  used  to  replace  the 
singularity-canceling  equation  at  maximum  load.  Such  a condi- 
tion, possibly  one  involving  the  amount  of  opening  at  the 
crack  tip,  must  therefore,  be  sought. 

(2)  At  the  minimum  load  level  in  a load  cycle,  the  nascent  super- 
dislocation is  decomposed  into  two  superdislocations  to  represent 
the  loss  of  plasticity  involved  in  unloading.  The  criterion 

for  the  decomposition  is  currently  one  that  maximizes  their 
strengths.  This  is  somewhat  arbitrary  and  certainly  other 
plausible  choices  exist  that  could  be  used.  The  computed 
growth  rates  are  sensitive  to  the  exact  choice,  however, 
particularly  when  K^n  = 0.  Thus,  a more  rigorous  procedure 
is  needed. 

(3)  The  super-superdislocation  concept,  in  which  the  superdisloca- 
tions generated  in  several  different  load  cycles  are  lumped 
together,  offers  both  significant  computational  efficiency  in 
the  complete  cycle-by-cycle  model  and  a simple  testing 
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ground  for  new  ideas.  Further  development  of  this  concept 
for  these  dual  purposes  is  therefore  warranted. 

To  take  the  model  as  it  presently  exists  to  the  point  where  the  objectives 
of  the  work  can  be  accomplished  will  require  improvements  in  the  manner  in 
which  the  computations  are  performed  as  outlined  above.  Measures  leading  to 
improved  qualitative  agreement  with  the  observed  R effect  and  crack-growth 
retardation  must  also  be  incorporated  into  the  model.  Comparison  with  alter- 
native approaches  and  with  crack  growth  and  other  experimental  data  (e.g., 
crack-closure  measurements)  can  then  be  made. 

To  improve  the  computational  efficiency  of  the  model,  a scheme  for  deter- 
mining when  and  how  the  superdislocations  representing  the  residual  plasticity 
from  previous  load  cycles  can  be  lumped  together  in  one  or  more  super-superdis - 
locations  should  be  pursued.  The  point  at  which  an  optimum  balance  between 
accuracy  and  computational  time  requirements  should  be  determined  (NB,  while 
this  point  could  depend  on  the  load  spectrum  under  consideration,  the  computed 
growth  rates  must  be  independent  of  the  manner  in  which  the  super-superdisloca- 
tions are  introduced.)  Linearization  of  the  governing  equations  to  achieve 
more  economical  computations  might  also  be  investigated. 

Finally,  it  must  be  emphasized  that,  while  the  work  performed  so  far  in 
developing  the  model  is  quite  promising,  there  can  be  no  assurance  that  the 
approach  described  above  will  ever  produce  a useful  tool  for  practical  applica- 
tions. It  is  possible  that  the  goal  of  a fundamentally  valid  model  that  can 
be  used  in  an  economical  manner  simply  cannot  be  attained.  If  the  inclined- 
strip  yield-superdislocation  model  cannot  be  constructed  to  meet  the  goal  of 
an  accurate  and  efficient  fatigue  crack-growth  predictive  technique,  then  such 
a goal  is  probably  not  within  the  current  state  of  the  art. 
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APPENDIX  A 


DERIVATION  OF  THE  BASIC  EQUATIONS  FOR  THE 
CYCLE-BY-CYCLE  ANALYSIS  MODEL  FOR  FATIGUE- 
CRACK  GROWTH  UNDER  ARBITRARY  SPECTRUM  LOADING 


General  Analytical  Approach 


The  analytical  approach  taken  in  this  work  rests  on  two  fundamental 
assumptions:  (1)  the  crack-tip  crack-opening  displacement  (CTCOD)  arising 

from  the  intense  plastic  deformation  attending  the  crack  tip  is  the  parameter 
that  governs  crack  growth  and  (2)  the  essential  feature  of  fatigue-crack 
growth  is  the  interaction  between  the  plastic  deformation  created  in  different 
load  cycles.  By  basing  the  work  upon  an  inclined  strip-yield  model,  the 
approach  taken  here  uses  what  is  perhaps  the  simplest  mathematical  model  in 
which  those  considerations  (i.e.,  CTCOD  and  distinguishable  plastic  zones  in 
each  load  cycle)  can  be  taken  into  account  in  a rigorous  fashion. 

Because  experimental  results  on  fatigue-crack  propagation  are  almost 
always  reported  in  terms  of  a linear-elastic  stress-intensity-factor  parameter, 
it  may  be  useful  to  show  the  connection  between  this  parameter  and  the  crack- 
opening displacement  (COD)  parameter  to  be  used  in  the  model.  Obtained  from 
the  Dugdale  strip-yielding  crack  model  [1]*,  the  commonly  used  expression  [2] 
for  the  latter  quantity 


- _ 8 YC  n a 

0 - - Y log  sec  2 Y 


(A-l) 


where  5 = COD  at  the  crack  tip 
2c  = crack  length 
a = applied  tensile  stress 
Y = yield  stress 
E = elastic  modulus. 

For  the  case  where  a is  small  in  comparison  to  Y,  the  above  expression  can  be 
linearized  to  get 

5 = . 

EY  " ' 


* References  for  this  appendix  are  given  on  page  78. 
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Now,  the  stress-intensity  factor  for  this  situation  is  simply  K = a/na . Hence, 
for  ct/Y  « 1,  it  can  easily  be  seen  that 


(A-2 ) 


This  means  that  for  the  load  levels  ordinarily  used  in  fatigue-crack-growth 
experiments,  treating  K as  the  driving  force  for  crack  extension  is  equivalent 
to  using  5 for  this  purpose.  Note  that  an  alternative  point  of  view  is  to  use 
the  Griffith  theory  to  express  the  strain-energy-release  rate  as  <?  = 6Y.  Then, 
from  the  equivalence  between  the  strain-energy-release  rate  and  stress-intensity 
factor  given  by  K2  = E£,  the  result  given  as  Equation  (A-2)  is  obtained. 

The  general  idea  upon  which  the  present  analysis  is  based--the  representa- 
tion of  a crack  and  its  attendant  zones  of  plastic  deformation  by  dislocation 
arrays --evolved  from  the  work  of  Head  and  Loust  [3]  and  Bilby,  et  al  [4]. 

When  both  approaches  are  tractable,  it  can  be  shown  that  the  dislocation  array 
approach  leads  to  solutions  that  are  identical  to  those  obtained  through 
continuum  elasticity,  e.g.,  the  Dugdale  problem.  The  virtue  of  this  dislocation 
approach  is  that  it  is  often  tractable  when  continuum  theory  is  not.  This  is 
true  in  particular  for  the  inclined  strip-yield  zone  solution  given  by  Bilby 
and  Swinden  [5],  which  is  the  basis  of  the  development  pursued  in  this  report. 

A drawback  of  the  dislocation  models  is  that  a singular  integral  equation  must 
be  solved.  This  would  make  the  Bilby-Swinden  model  per  se  exceedingly  cumber- 
some for  use  in  representing  fatigue-crack  growth  despite  its  advantageous 
features.  Consequently,  the  simplification  afforded  by  Atkinson's  superdislo- 
cation concept  [6-8]  has  been  introduced  to  enable  a large  number  of  load 
cycles  to  be  examined. 

The  notation  used  in  the  following  generally  corresponds  to  that  used  by 
Rice  [9],  Results  not  specifically  referenced  elsewhere  can  be  found  in 
Reference  9. 

Mathematical  Development 

It  will  be  considered  in  this  report  that  crack  growth  occurs  in  an 
isotropic  linear  elastic  body  under  plane  strain  (or  generalized  plane  stress) 
conditions.  Crack-tip  plasticity  will  be  represented  by  discrete  dislocations, 
finite  in  number,  that  are  confined  to  specified  slip  planes.  Because  discrete 
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dislocations  can  be  treated  within  the  theory  of  elasticity,  it  is  appropriate 
to  adopt  the  complex  variable  approach  developed  by  Muskhelishvili.  In  this 
approach,  the  stresses  and  displacements  for  elastic  deformation  in  the  xy 
plane  are  expressed  in  terms  of  two  arbitrary  functions  of  the  complex 
variable  z = x + iy  bv  the  relations 


Here,  G is  the  shear  modulus  and 


where  v is  Poisson's  ratio.  Note  that  in  Equations  (A-3)  through  (A-5)  and  in 
the  following,  the  prime  notation  indicates  differentiation  with  respect  to  z, 
the  bar  notation  is  used  to  denote  the  complex  conjugate,  and  the  convention 
F(z)  = F(z)  is  used  throughout.  Also  note  that  + and  - signs  used  as  subscripts 
denote  limits  of  the  function  as  the  x-axis  is  approached  from  the  regions 
y > 0 and  y < 0,  respectively. 

The  basic  relations  needed  are  those  for  the  interaction  between  a dislo- 
cation and  a crack  in  an  infinite  body  under  remotely  applied  loading.  It  is 
considered  that  the  crack  lies  along  the  x-axis  with  the  tips  of  the  crack 
being  at  the  position  (a,0)  and  (-a,0).  The  dislocation  is  located  at  the 
position  Z = § + iT]  and  has  a Burger's  vector  B = bx  + ib2 . Although  we  will 
confine  our  attention  to  tensile  loading  in  treating  crack  growth,  at  this 
point  in  the  development  of  the  model  a more  general  applied  loading  of  biaxial 
tension  plus  shear  is  considered;  i.e.,  as  f xa  + y3  - ® c?11  = xx  , a22  = T2,  and 
CT12  = S (see  Figure  A-l). 

A solution  for  the  stress  and  displacement  fields  for  the  problem  shown 
in  Figure  A-l  can  be  obtained  by  linear  superposition  after  first  determining 
the  complex  potentials  for  two  subproblems.  These  subproblems  correspond  to 
an  isolated  dislocation  and  an  isolated  crack. 


FIGURE  A- I.  LOADING  CONFIGURATION 


The  Isolated  Dislocation 


The  complex  potentials  for  a dislocation  with  Burger's  vector  bx  + ib 
at  the  generic  position  Z = § + iT|  are  given  by 

G (b1  + ib2 ) i 

/ / \ 1 ^ / A "7  \ 


where 


Combining  Equations  (A-3)  through  (A-5)  with  (A-9)  gives 


Or,  by  substituting  Equations  (A-7)  and  (A-8)  into  the  above,  we  have 


for  the  dislocation. 

It  is  important  to  recognize  that  here  and  throughout  this  work,  care 
must  be  taken  with  the  branches  of  the  log  function  in  the  expression  for  the 
displacements.  It  is  this  function,  of  course,  that  gives  the  discontinuitv  in 
displacement  characteristic  of  the  dislocation.  Relations  such  as  Equation 
(A-12),  it  might  be  noted,  can  be  used  in  incorporating  the  crack-closure 
effect  [10-12]. 


For  the  crack  problem,  it  can  be  shown  that  for  the  stresses  to  vanish 
at  infinity,  we  must  have 


Equations  (A-3)  through  (A-5)  can  then  be  written  in  terms  of  a single  potential 
function  as  follows: 


CT11  + a22  = 4 Re  [a  '(z)}  > 

k (oa2  -all>  + ia!2  = a'(z)  - - (z  - z)  ®"(z)  , 


(A- 14) 


2G  (uT  + iu2)  = Ka(z)  - (3(z)  + (z  - z)  £ (z) 


(A- 15) 


For  these  conditions,  Rice  shows  that  the  potential  function  is  given  by  the 


relation 


*'(*>  - - ^ (**  -azyk  J tPg (t)  - ipx  (t)  ] fe(V-z)?%  dt  • 


(A-16) 


The  simplest  situation  is  where  the  surface  tractions  are  constant.  This 
arises,  in  particular,  when  the  solution  for  a crack  in  a body  acted  upon  by 
uniform  remote  tension  and  shear  is  to  be  obtained  by  superposition.  Referring 
to  Figure  A-l,  the  applied  load  exerts  a stress  field  on  the  x-axis  that  is 

ct22  “ a = - iS 

Hence,  we  have  p,  = - S2  and  p2  = T2.  Substituting  these  into  Equation  (A-16) 
then  gives 


Ta-  IS 


so  that 


(z)  = - \ (T2  - i3)  1 


(z3  - a3)' 


It  is  easily  found  that 


a"  (z)  = - k (t,  - is)  ■ -q  2 , 

(z3  -a3) j 


and,  by  integrating  Equation  (A-17), 


(z)  = - k (t2  - is)  [z  - (z3  -a3)^]  . 


(A-17) 


(A- 18) 


(A-19) 


The  constant  of  integration  corresponds  to  rigid  body  translation  and  has  been 
omitted  from  Equation  (A-19). 

Because  (z~  -2e)^  = i (2‘  - z-1)^  and  taking  the  positive  value  of 
(a3  -z3)*  on  the  upper  side  of  the  branch  cut,  it  is  found  that  when  y - 0, 
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shkd* 


= (iT2  +s)  (a“  -iry  , 

on  the  interval  |x|  < 2 and  zero  outside  the  interval.  Hence,  from  Equation 
(A-15),  the  discontinuity  in  the  displacements  on  the  crack  faces  due  to  the 
applied  stress  can  be  obtained  from 


[ux  + iu2l  - [uj  +iu2]_  - 


(iT,  +S)  (a3  -xe) 


Hence,  the  crack-opening  displacement  is  given  more  explicitly  by 
Aiii  + iAu2  = (S  + iT-,)  - x2 ) 2 |x|  S3  , 


(A-20) 


To  obtain  the  displacement  elsewhere  in  the  body,  the  uniform  field  due  to 
stresses  cl:  = Tx  , a22  = T2  and  cr12  = S must  be  superposed.  Using  Equations 
(A-3)  through  (A-5),  this  is  easily  found  to  give 

u,  + iu2  = ^ (Tx  +T2)  z+\  (Tx  - T_  + 2iS ) zj  . 

But,  because  this  expression  gives  zero  contribution  to  the  COD  under  tensile 
loading  alone,  it  will  not  be  used  in  the  following. 


Dislocation  and  Crack  Together 


To  obtain  a solution  for  a dislocation  located  near  a crack,  we  first 
determine  the  stresses  exerted  by  an  isolated  dislocation  on  the  line  segment 
where  the  crack  is  to  be  located.  For  the  dislocation  positioned  as  shown  in 
Figure  A-l,  Equations  (A-10)  and  (A-ll)  should  be  combined  to  get  an  expression 
for  the  normal  and  shear  stress  exerted  on  the  x-axis.  This  is 


r G(bx+ba)  \ G r 2bx 

+ ict12  = 2Re  ;iTT  (R  + 1}  (x  . ?)/  'in  (X  + 1)  '„x  - Z 


+ (A_n) 


(x  -7Y 


The  next  step  is  to  determine  the  potential  function  for  a crack  with  stresses 
acting  on  its  surfaces  which  will  wipe  out  the  dislocation  stresses  on  the 
crack.  This  is  accomplished  by  substituting  Equation  (A-21)  into  Equation  (A- 
16).  To  integrate  the  result,  we  need  the  following  general  result  which  can 
be  obtained  by  contour  integration: 


(c  - z)  (t  - Z) 


z - Z 


As  above,  the  branch  of  the  complex  square  root  is  defined  such  that  (aa  - z 
is  the  branch  cut  along  (-a,  a)  taking  the  positive  value  on  the  y > 0 side 
Also,  note  that  for  large  z,  (ac  -z2)-  - - iz  + 0(1)  so  that  (z2  -a3 *)'2  = 
i (CZ2  -z2)^.  With  these  results  we  obtain 


2ni  (K  + 1)  (zs  _ z2  ) 


(Z*  (Z  - z)2 


It  will  be  necessary  in  the  following  to  have  an  explicit  expression  for 
(z)  as  well  as  (z) . This  is  obtained  by  differentiating  Equation  (A-23). 
The  result  is 

„ G r -z  r (Z5  (Z2  -a2)^- 

d (Z)  = 2ni  (K  -Hi)  (bi+lbc)  l(z2  -a2)a  -2  + z - Z + z - Z _ 


1 r(z 


(za  -a2)' 


2 - -a  ) ^ - (Z 
(z  - Z) 


:2  -s2)^  (z2  - -2  r?  - (z2  - -2 ) -~ 

2 (z  - Z)2  - J 


(A-24) 


G (Z  - Z) (Z  - z ' 

2ni  (K+l)  1 12  (zs  -s2)^  (Z  -C2)^  (Z  - z)2  lz3  - 22 


2 (zB  -2s  )%  (Zc-32)^  + 2s2  - 2Zz 

Z - z / 

This  completes  the  derivation  of  the  basic  equations.  We  now  turn  to  the  inclined 
strip-yield  model  itself. 

The  Inclined  Strip-Yield  Crack  Model 

Basis  for  the  Model 

The  equations  developed  in  the  preceding  section  describe  the  interaction 
between  a crack  and  a dislocation  in  an  infinite  isotropic  body  with  remotely 
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applied  loads  under  plane-strain  conditions.  These  are  the  fundamental  ele- 
ments to  be  used  in  formulating  the  solution  for  a crack  with  crack-tip  plas- 
ticity represented  by  strip-yield  zones.  As  shown  in  Figure  A-2,  we  intend  to 
characterize  the  extent  of  yielding  by  the  position  and  strength  of  a single 
superdislocation  on  each  slip  plane.  Note  that  the  angle  between  the  slip 
plane  and  the  crack  plane  is  not  to  be  associated  with  a particular  crystallog- 
raphic direction.  Instead,  the  intent  is  to  characterize  plastic  deformation 
over,  perhaps,  several  grains.  This  means  that  the  angle  should  be  at  least 
45  degrees  to  be  in  general  accord  with  the  shape  of  macroscopic  plastic-yield 
zones  observed  around  crack  tips;  see,  for  example,  the  experimental  evidence 
presented  by  Hahn,  et  al  [13]. 

In  this  section,  the  development  is  confined  to  a single  set  of  slip 
planes  as  would  occur  in  one-half  cycle  of  loading.  It  is  intended  that  the 
deformation  be  symmetric  with  respect  to  the  origin  of  coordinates.  Conse- 
quently, it  will  suffice  to  describe  the  location  and  strength  of  the  disloca- 
tion in  the  first  quadrant  (i.e.,  position  1 in  Figure  A-2).  As  shown  in 
Figure  A-3,  we  will  let  2 a be  the  crack  length,  Z be  the  distance  from  the 
crack  tip  to  the  dislocation,  8 be  the  angle  between  the  crack  plane  and  the 
slip  plane,  and  b be  the  strength  of  the  dislocation. 

If  Bj  denotes  the  Burger's  vector  of  the  dislocation  at  position  Z , it 
will  be  convenient  to  let 

Bj  « bji  , J « 1,  2,  3,  4 , (A-25) 

where  the  jj.  1 s are  dimensionless  quantities  which  depend  only  on  0.  We  further 
note  that  the  positions  of  all  of  the  dislocations  can  be  expressed  relative  to 
the  position 

Z1  = a + 2e10  . (A-26) 

For  references  purposes  in  the  following  discussion,  it  is  convenient  to  summa- 
rize the  locations  and  strengths  of  the  dislocations  as  in  Table  A-l. 


TABLE  A-l.  LOCATIONS  AND  STRENGTHS  OF  DISLOCATIONS 
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The  stress  and  displacement  fields  for  the  inclined  strip-yield  problem  are 
obtained  by  superposition  using  the  solutions  for  the  following  subproblems: 

• Four  symmetrically  placed  isolated  superdislocations 

• A crack  with  nonuniform  surface  tractions  equal  and  opposite  to 
the  dislocation  stresses 

• Remotely  applied  tension 

• Crack  with  uniform  traction  equal  and  opposite  to  the  applied 
stress. 

In  each  of  the  two  above  crack  solutions,  there  is  a singularity  at  the  crack 
tips.  It  is  taken  as  axiomatic  in  this  work  that  the  combined  singularity  must 
be  abolished.  This  proviso  gives  one  equation  which  can  be  used  to  help  deter- 
mine the  two  unknowns  occurring  in  the  problem.  The  two  unknowns  are  the 
strengths  and  positions  of  the  dislocations;  i.e.,  the  quantities  b and  i.  The 
second  equation  needed  is  obtained  by  assuming  that  the  net  resolved  shear 
force  acting  on  a dislocation  in  the  slip  plane  is  balanced  by  a known  retarding 
or  friction  force  related  to  the  yield  stress. 


The  Singularity-Cancelling  Equation 

The  singularity  at  the  crack  tips  can  be  abolished  by  setting  the 
coefficient  of  the  singular  term  in  the  stress  field  equal  to  zero.  More 
directly,  by  combining  Equations  (A-13)  and  (A- 14),  we  find  that 

ct22  * i-CTi  2 = $ ' (z)  + 3 (z)  + (z  -z)  <&"  (z)  . 

Consequently,  on  the  x-axis, 


t . • 


cr22  - ia12  = t+'  (x)  + tj  (x)  . 

Or,  using  Equation  (A-17)  with  T2  = cr  and  S = 0 together  with  Equation  (A-23), 
we  obtain 


*-'22  ” : 


a h -■ 


(X2  -a3)^  TT  (K  + l)  '-j-l  (X2  -23)$ 

(z  s -s2)^  - (x2  -a2)^  (Z.2  -a2)^  - (x2  -s2)3?. 


rtf  -a* 


(Z;2  -22)^+32  - Z ,x- 


(z.  - xf  (!  * -22)^ 


X . 


t *-r  *> 

nfVtMrtr  ■ " i-ii 


J ’ll.": 


The  singularity-cancelling  equation  can  be  extracted  from  this 
expression  almost  by  inspection.  It  is  the  coefficient  of  the  term  (xc  -2s)' 
evaluated  at  x = 2 . By  setting  this  to  zero,  we  obtain 


0 =cfl  - • 


(K  + l)  l J 


,3  \ 2 / r,  3_~2\  2 


{z.a -a  ) (z.-2)i  4 

+— i — = — + ■ 1 ^ l+2  s 


-<*>n  4 ^ d - Z ) 

I CL  Y.  ■■  ■ ' J — J ■ i 

J j=i  a.  -2)  (z 2 - 2s ) ; 

J J j j 


From  the  data  given  in  Table  A-l,  it  is  evident  that 


T.  Li, . =0 

j=l  J 


Consequently,  the  singularity-cancelling  equation  reduces  to  an  expression  that 
can  be  written 

4 (Z  3 - )%  f~Z  3 -'z3)% 

Gb  r 2 i “ ’ ' i 1 ' 

O'  2 = — — , , , Im  1.  u , — f + — = 

rr  (K  + l)  C . , j _ Z.-2  Z.-2  J 

J=1  J J 


(A-27 ) 


(Zi  - ZJ> 


+ 4 E [l.  — 4 _ J 1 

j=l  J (zj  -3)  (z. -32)^  ; 


Equation  (A-27)  is  the  first  of  the  two  basic  equations  needed  for  determining 
the  two  unknowns  b and  l.  Notice  that  this  relation  is  to  be  interpreted  for 
given  values  of  the  applied  stress  and  crack  length,  for  a single  application 
of  the  load  only. 
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APPENDIX  B 


EXPERIMENTAL  MEASUREMENTS  OF  FATIGUE-CRACK 
GROWTH  AND  RETARDATION  IN  Ti-6AI-4V  ALLOY 


Introduc Cion 


Selected  fatigue-crack-propagation  experiments  were  conducted  to  provide 
substantiation  of  the  analytical  developments  of  this  research  program.  Included 
were  the  generation  of  fatigue-crack-propagation  data  under  conditions  of 
constant  cyclic  stress-intensity-factor  range  (i.e.,  constant  AK) , the  study  of 
the  overload  effects  on  crack-growth  retardation,  the  calibration  of  specimen 
compliance,  and  observations  of  crack-closure  effects.  This  experimental  effort 
was  directed  toward  complementing  and  refining  results  obtained  on  previous 
NADC-sponsored  programs  and  to  establish  a firmer  experimental  basis  for  the 
iterative  approach  to  this  experimental-analytical  development  activity. 

Experimental  Approach 

Fatigue-crack-propagation  data  are  usually  generated  using  test  specimens 
which  directly  simulate  cracked  structural  elements  or  which  satisfy  the  con- 
straints of  fracture-toughness  testing.  Under  constant -amplitude  loading, 
these  specimens  exhibit  increasing  crack-growth  rates  with  advancing  crack 
length  since  the  applied  stress-intensity  factor  is  also  increasing.  However, 
in  this  program,  to  achieve  a more  definitive  description  of  crack  behavior, 
fatigue-crack-propagation  data  were  generated  under  conditions  of  constant 
cyclic  stress-intensity-factor  range,  AK.  In  addition  to  providing  the  basic 
data  for  quantifying  the  free  parameters  of  the  evolving  model,  this  approach 
provides  a clearer  picture  of  the  influence  of  the  dominant  mechanical  variables 
for  several  reasons.  First,  if  the  rate,  da/dN,  is  uniquely  dependent  upon  the 
stress-intensity  factor,  K,  or  on  its  range,  AK,  the  crack-growth  curve  of 
crack  length  versus  cycle  count  will  be  linear.  This  would  permit  a direct  and 
convenient  means  of  fitting  crack  growth  which  is  complexly  nonlinear  under 
the  more  usual  conditions  of  constant-amplitude  loading.  Second,  the  dependency 
of  fatigue-crack-propagation  rate,  dA/dN,  on  K or  AK  can  be  more  clearly 


established.  Third,  a large  number  of  overload  combinations  can  be  studied 
over  a range  of  crack  lengths  in  one  specimen. 

Although  a constant  cyclic  stress-intensity-factor  range,  AK,  also  can  be 
achieved  through  configuration  control,  as  in  the  tapered  double-cantilever 
beam  (DCB)  specimen,  such  geometric  constraints  are  not  generally  desirable. 
The  tapered  DCB  specimen,  while  available  to  provide  a constant  AK  independent 
of  crack  length  under  constant-amplitude  loading,  does  not  simulate  the  broad 
class  of  crack  configurations  encountered  in  structures  nor  is  there  a large 
body  of  such  fatigue-crack-propagation  data  for  correlation  of  the  evolving 
model. 


Program  Details 


The  basic  experimental  details  of  this  research  program  are  presented  in 
the  following  subsections.  Descriptions  of  materials,  test  specimens,  and 
experimental  procedures  are  included. 


Materials 

The  fatigue-crack-growth  and  retardation  studies  were  conducted  on  both 
mill-annealed  and  beta-processed  Ti-6A1-4V  titanium  alloy  plate.  The  mill- 
annealed  material  (normally  ^-inch  thick)  produced  by  Reactive  Metals,  Inc. 

(RMI),  was  available  at  Battelle  from  a previous  NADC- sponsored  research  pro- 
gram. The  beta-processed  material  (normally  5/16-inch  thick)  produced  by 
Titanium  Metals  Corporation  of  America  (TMCA)  was  provided  by  NADC  from  surpluses 
in  the  Federal  Aviation  Administration's  Supersonic  Transport  Program  at  the 
Boeing  Company. 


Chemistry.  The  average  chemistry  for  these  materials  is  presented  in 
Table  B-l.  The  mill-annealed  plate  chemistry  was  certified  by  RMI  for  Ingot 
295338.  The  beta-processed  plate  chemistry  was  certified  by  TMCA  for  Heat  No. 
K2701  and  verified  by  Boeing  Company. 
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TABLE  B-l.  CHEMICAL  COMPOSITION 


Weight,  percent 


Material 

A1 

V 

Fe 

0 

C 

N 

H 

Mill  Annealed 

6.4 

4.2 

0.18 

0.127 

0.022 

0.010 

0.003 

Beta  Processed 

6.4 

4.2 

0.07 

0. 150 

0.022 

0.028 

0.0060 

Mechanical  Properties.  The  static-tensile  properties  of  each  material  are 
presented  in  Table  B-2.  The  properties  of  the  mill-annealed  plate  were  deter- 
mined at  BCL;  those  of  the  beta-processed  material  were  developed  by  TMCA  and 
verified  by  Boeing. 


TABLE  B-2.  MECHANICAL  PROPERTIES 


Mill  Annealed 
Beta  Processed 
Beta  Processed 


L 

137 

130 

15 

35 

L 

139 

126 

17 

29 

r 

140 

126 

17 

30 

Grain 

Direc- 


Tensile 
Ultimate 
Strength, 
TUS,  ksi 


0.2  Percent 
Offset  Tensile 
Yield  Strength, 
TYS,  ksi 


Reduction 
in  2-Inch  in  Area, 
Gage  Length,  RA, 
e,  percent  percent 


Metallography.  The  metallography  of  the  mill-annealed-plate  material  has 
been  discussed  in  the  final  report  of  the  NADC-sponsored  research  program.  The 
following  discussion  describes  the  metallographic  characterization  of  the  beta- 
processed  material. 

Sampling  of  the  beta-processed  Ti-6A1-4V  plate  was  accomplished  to  enable 
examination  in  three  planes--surface,  longitudinal  thickness  (parallel  to  roll- 
ing direction),  and  transverse  thickness  (transverse  to  rolling  direction). 
Standard  metallographic  techniques  were  used  to  grind  and  polish  samples. 
Krolls'  etchant  (HF,  HN03 ) was  used  to  reveal  microstructure.  Samples  were 
reetched  to  show  more  clearly  the  grain  size  variation  that  exists  across  the 
plate  thickness.  Representative  photomicrographs  are  shown  in  Figures  B-l  and 
B-2. 
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FIGURE  B- 1 . THICKNESS  CROSS-SECTION  OF  BETA-PROCESSED  TITANIUM  PLATE 
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FIGURE  B-2.  THREE-DIMENSIONAL  VIEW  OF  MICROSTRUCTURE 


The  large  grain  size  of  the  prior  beta  grains  in  this  plate  of  transformed 
beta  microstructure  is  a feature  almost  discernible  with  the  naked  eye.  As 
shown  in  Figure  B-l(a),  the  cross-section  photomacrograph  reveals  another 
feature--the  variation  in  prior  beta  grain  size  that  exists  between  the  more 
heavily  worked  surface  layers  and  center  material.  The  smaller  surface  layer 
grains  are  about  300  microns  across  (ASTM  Grain  Size  about  1),  while  the  grains 
of  the  center  material  are  three  to  eight  times  larger. 

Grain  size  variation  across  plate  thickness  is  not  an  unusual  feature  in 
titanium  alloys  but  as  shown  in  Figure  B-l(a),  the  structure  of  this  plate 
is  characterized  by  an  unusually  abrupt  change  in  grain  size  between  surface 
layers  and  center  material.  As  shown  in  Figure  B-l(b)  and  (c),  this  is  not  a 
large  difference  between  the  microstructures  of  the  two  different  grain  size 
materials.  The  acicularity  of  the  alpha  platelets  is  a bit  more  defined  in 
the  center  material  but,  basically,  the  transformation  products  and  form  are 
the  same  from  both  small  and  large  grains.  A feature  of  both  microstructures 
is  the  lack  of  a pronounced  alpha  lining  at  prior  beta  grain  boundaries.  Where 
alpha  grain  boundary  lining  is  apparent,  it  is  thin,  discontinuous,  and  not  well 
developed,  indicating  the  preferred  degree  of  rapid  cooling  from  the  beta 
temperature. 

A pronounced  elongation  of  grains  in  the  rolling  direction  is  not  observed 
in  this  material  although  a slight  alignment  of  the  prior  beta  grains  is  apparent 
(Figure  B-l).  This  lack  of  fibering  is  revealed  further  in  the  three-dimensional 
view  of  the  microstructure  as  presented  in  Figure  B-2  (100X).  Again,  the  out- 
standing feature  of  this  transformed  beta  structure  is  the  prior  beta  grain  size 
difference  between  surface  and  center  material.  Differences  in  mechanical 
properties  might  be  expected  between  materials  (or  across  thickness  sections) 
with  this  much  difference  in  grain  size. 

Test  Specimens 


The  specimen  design  used  in  the  program  is  illustrated  in  Figure  B-3.  The 
single  large  pin  hole  drilled  at  each  end  was  used  to  provide  more  positive 
alignment  in  these  relatively  low  load  tests.  The  actual  i-lnch  hole  was  saw- 
notched  as  a crack  starter. 


Experimental  Procedures 

All  tests  were  conducted  in  BCL's  Structures  Laboratory  on  a 50-kip  servo- 
controlled  electrohydraulic  testing  system. 

Fatigue-Crack  Propagation.  The  central  hole  of  each  test  specimen  was 
notched  with  a fine  saw  cut  and  precracked  to  an  initial  crack  length  of  approx- 
imately 0.50  inch.  Each  specimen  was  then  cyclically  loaded  over  a stress  range 
to  provide  the  desired  AK  and  R values.  Surface  crack  length  was  monitored  with 
an  optional  comparator.  Crack  length  measurements  were  made  at  approximately 
every  0.050  inch  at  which  time  the  loading  was  decremented  to  maintain  the 

desired  K value  according  to  the  re'ition 
max 

a = K [W  tan  -2/W]‘‘  , (B-l) 

where  a = gross  applied  stress 

K = applied  stress-intensity  factor 
W = panel  width 
a = half-crack  length. 

The  crack  length  and  associated  cycle  counts  were  recorded  and  plotted  to  main- 
tain a continual  record  of  crack  growth. 

In  summary,  the  transformed  beta  microstructure  of  this  plate  appears  to 
be  typical  of  Ti-6A1-4V  betatized  material.  The  differences  between  the  prior 
beta  grain  sizes  of  surface  and  center  areas,  as  well  as  the  abrupt  change  in 
grain  size  between  areas,  is  a noteworthy  feature. 

Overload  Procedures.  The  overload-crack-retardation  experiments  were 
conducted  as  a simple  extension  of  the  procedures  described  above.  Single  over- 
loads (or  overstresses)  were  applied  periodically  during  the  basic  fatigue- 
crack-propagation  tests.  A load-compliance  record  was  made  for  each  overload 
such  that  increases  in  crack  extension  could  be  related  to  specimen  compliance. 

A miniaturized  clip  gage  was  mounted  in  the  mouth  of  the  EDM  notch  to  monitor 
specimen  compliance  during  loading. 

Matrix  of  Test  Data 

Crack  propagation  in  these  two  materials  was  evaluated  both  with  and  with- 
out overloads.  The  matrix  of  data  generated  is  indicated  as  follows: 


Maximum  Cyclic 
Stress-Intensity 
Factor,  K^, 
ksi-in.2 

Stress 

Ratio, 

R 

Number  of  Test  Specimens 

Mill 

Annealed 

Beta  Annealed 

Without 

Overload 

With 

Overload 

Without 

Overload 

With 

Overload 

10 

0.05 

1 

1 

-- 

-- 

20 

0.05 

2 

1 

1 

1 

20 

0.05 

2 

1 

1 

1 

number  of  overloads 

and  overload  ratios  were  applied 

to  each  of 

the  overload 

specimens.  Retardation  was  considered  to  have  been  overcome  when  the  crack- 
propagation  rate  as  determined  under  constant  AK  conditions  was  reestablished. 


Experimental  Results 

The  principal  experimental  results  derived  in  this  research  program  are 
summarized  in  the  following  sections.  While  these  data  provide  an  initial  basis 
for  the  analytical-experimental  iteration,  they  will  be  expanded  as  this  model 
development  is  continued. 


Constant  AK  Crack-Growth  Behavior 


A typical  crack-growth  curve  generated  under  conditions  of  constant  s tress  - 
intensity-factor  range  is  shown  in  Figure  B-4.  The  curve  has  the  expected 
linearity  until  the  crack  length  exceeds  half  the  panel  width  whereat  some 
divergence  is  noted.  At  maximum,  the  curve  slope  prior  to  fracture  diverged 
from  the  initially  established  curve  by  a factor  of  3. 

The  cause  of  this  divergence  could  be  attributed  either  to  an  inadequacy 
of  the  width  correction  in  defining  the  AK  level  of  loading  or  to  other 
mechanical  factors  affecting  the  growth  rate.  In  the  light  of  the  extensive 
mechanics  analyses  on  this  specimen  configuration,  the  former  case  is  highly 
unlikely.  Therefore,  it  is  considered  that  other  net  section  stress  influences 
come  into  play. 

Fatigue-Crack-Propagation  Rate 


The  fatigue-crack-propagation-rate  behavior  derived  from  these  conducted 
under  conditions  of  constant  cyclic  stress-intensity-factor  range  are  summarized 
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IGURE  15-4.  TYPICAL  CONSTANT  AK  CRACK-GROWTH  CURVE  (SPECIMEN 


h 


in  Figure  B-5.  The  results  are  presented  over  dual  abscissa  formats.  On  the 
left,  the  rates  are  presented  in  terms  of  the  actual  stress-intensity-factor 


range, 


AK  = (1  - R)  K 


(B-2) 


on  the  right,  the  results  are  presented  in  terms  of  effective  stress-intensity- 


factor  range, 


AK  ,,  = (1  - R)  K 
e f f s max 


(B-3) 


where  m is  the  Walker  [1]*  exponential  parameter  derived  from  a collapse  of 
test  data  on  the  mill-annealed  material.  This  display  shows  both  the  effect  on 
stress  ratio  and  the  difference  in  behavior  of  the  two  materials  tested. 

It  should  be  noted  that  these  rate  results  are  presented  in  terms  of  total 
crack  length,  2a,  for  consistency  in  comparison  with  previously  published 
results.  In  the  main  body  of  this  report,  the  analytical-experimental  compari- 
sons have  been  presented  in  terms  of  the  half-crack  length,  a,  for  consistency 
with  the  generalized  model, 

A broader  perspective  of  these  results  may  be  gained  by  comparing  them 
with  the  summary  results  of  a previous  NADC-sponsored  program  [2]  as  shown  in 
Figure  B-6.  The  shaded  areas  represent  previous  test  data  generated  under  condi- 
tions of  constant-amplitude  loading.  The  solid  lines  represent  the  average  rate 
results  of  these  data  modeled  by  the  formulation 


d(2a)  _ 


(1  - R)  Kc  - AK 


(B-4) 


M 

I.  J-.: 


where  C = 0.267  u,in. /cycle 
n = 3.30 

t 

K = 267  ksi-in. 
c 

The  results  of  this  present  program,  translated  directly  from  Figure  B-5,  are 
overlayed  on  the  previous  results  as  individual  symbols.  The  correlation 
appears  to  be  quite  good. 

A difference  in  crack  behavior  in  the  two  materials  was  also  notable  in 
the  macroscopic  appearances  of  their  crack  surfaces.  As  shown  in  Figure  B-7, 


* References  for  this  appendix  are  listed  on  page  100. 


i a 

l ' 


V*T  ’?-• 


W 


Propagation  Rate,  d(2a)/dN,  inch/cycle 


Material 


Beta  Mill  Stress 
Processed  Annealed  Ratio,  R 


FIGURE  B-5 . CRACK-GROWTH- RATE  BEHAVIOR 
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FIGURE  B-7.  CONTRASTING  APPRARANCES  OF  CRACK  SURFACES 


the  BCL  material  which  exhibited  the  higher  rate  behavior  had  a relatively 
smooth  "satiny"  finish;  whereas,  the  NADC  material  which  exhibited  a lower  rate 
had  a coarser  textured  finish.  This  difference  is  attributed  to  the  beta 
processing  of  this  material. 

Influence  of  Overloads  on  Retardation 

Retardation  was  measured  as  the  additional  number  of  cycles  required  for 
a crack  to  grow  from  the  point  of  overload  to  a recovery  point  at  which  the 
previously  monitored  uninterrupted  crack-propagation  rate  was  reestablished. 

This  is  determined  from  the  basic  crack-growth  curve  as  the  horizontal  displace- 
ment of  tangents  of  equal  slope  before  and  after  the  overload  application. 

A summary  of  the  overload  test  results  is  presented  in  Figure  B-8.  The 
number  of  kilocycles  of  retarded  crack  growth,  NR,  is  plotted  as  a function  of 
the  percent  of  single  overload  applied.  Although  there  is  considerable  scatter 
in  the  data,  a distinct  influence  of  stress  ratio  is  evident.  The  number  of 
cycles  retarded  increases  with  stress  ratio,  but  appears  to  be  independent  of 

the  K value  over  the  range  studied.  The  retardation  behavior  appears  to  be 
max 

equivalent  for  both  materials,  although  the  general  crack-propagation  rates  were 
different. 

Compliance  Calibration 


To  demonstrate  and  verify  the  elastic  response  of  the  test  specimens,  speci- 
men compliance  was  calibrated  against  crack  length.  The  clip  gage  located  in 
the  central  hole  was  used  to  monitor  displacement  adjacent  to  the  crack  surface. 
Load  deflection  records  were  made  after  precracking,  during  the  overload  applica- 
tions, and  at  other  crack  lengths  during  the  crack-growth  tests.  The  following 
analytical  formulation  of  specimen  compliance  was  used. 

For  a system  of  colinear  cracks  under  a remote  tensile  stress  field,  the 
Westergaard  [3]  stress  function,  Z,  and  its  first  integral,  Z,  may  be  expressed, 
respectively,  as 


Z 


a ri 


■■na/W\2 

CT  J 


( B - 5 ) 


and 
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Percent  Overload 


- jW  _i  /cos  nz/W\  . _ . aW  , _i  cos  ttz/W\  . _ 

Z = — cos  1 —777  + C = l — cosh  1 —777 ) + C . (B-6) 

tt  vcos  rra/Wy  rr  \cos  na/W/  v 

Since  the  relative  crack-surface  displacement,  2T|0  , in  the  y direction  may  be 
expressed  as 


2T|0  = 


4(1  - y3 


Im  Z , 


it  follows  that 


2T1  = .^0_lAL3 ) SZ  C0sh-1  'cos  mc/W\ 
iT1°  E tt  costl  Vcos  m/w/ 


(B— 7 ) 


(B-8) 


for  -a  £ x 5 a.  For  the  paeticular  case  of  x = 0,  the  centerline  crack  surface 
displacement,  A,  may  be  expressed  as 


2Tlo  , x = 0 = A = 


4(1  - v3)  <jW 


aw  ,_i  ( rns\ 

— cosh  ; sec  — 


( B— 9 ) 


Accordingly,  the  stress-normalized  effective  compliance  may  be  expressed  as 


A 4(1  - v2 


W cosh-1  (sec  rra/W) 


(B- 10) 


A comparison  of  analytical  and  experimental  results  based  on  a panel  width 
of  W = 6 inches  and  a material  modulus  of  E = 15  x 10s  psi  is  presented  in  Figure 
B-9.  The  correlation  appears  very  good  and  confirms  the  adequacy  of  the  analytical 
formulation. 

Observations  on  Crack  Closure 


Fatigue-crack  closure,  a phenomenon  first  documented  by  Elber  [4],  is  the 
closure  or  impingement  of  fatigue-crack  surfaces  prior  to  the  complete  unloading 
of  a test  specimen  containing  a fatigue  crack.  If  one  considers  the  load-deforma- 
tion record  of  a test  specimen  wherein  the  deformation  is  monitored  in  the  crack- 
tip  region,  an  initial  nonlinearity  will  be  observed  although  the  overall  loading 
and  structural  conditions  are  linearly  elastic.  This  is  illustrated  in  Figure 
B-10(a).  The  nonlinearity  is  attributed  to  the  compression  (or  relaxation)  of 
the  crack  surfaces  as  closure  occurs  (or  is  relieved)  during  the  loading  cycle. 

From  this  figure  it  can  be  seen  that  the  effectivity  of  closure  is  dependent 

upon  its  relationship  to  the  maximum  and  minimum  cyclic  stresses  respectively, 

a and  a . . That  is,  if  the  opening  stress,  a.,  is  above  a . , the  crack 
max  min  * r ’o’  min’ 
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FIGURE  B-9.  COMPLIANCE  CALIBRATION 
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CRACK-CLOSURE  BEHAVIOR 


1 


will  not  be  opened  in  full  proportion  to  the  applied  stress  range,  a ” amin' 

Rather  its  opening  and  associated  damage  will  be  attenuated  in  proportion  to 

a - a . This  gives  rise  to  the  effective  stress-range  ratio  defined  as 
max  o ° ° 


U 


a 

max 


a 

max 


for  a0 


> a . 
min 


(B-ll) 


where  amax  and  are  the  stress  extremes  of  the  preceding  fatigue  cycles. 

Where  U is  greater  than  unity,  closure  is  not  considered  to  influence  crack 

growth.  Where  U is  less  than  unity,  closure  is  considered  to  be  a direct  scalar 

factor  on  AK  in  crack-growth  calculation. 

Crack-opening  stresses,  aQ  , were  determined  by  the  use  of  an  Elber-type 

closure  gage  at  the  crack  tip  as  well  as  from  the  records  used  for  calibration 

on  specimens  without  overloads.  The  resultant  values  are  generally  normalized 

in  Equation  (B-ll).  The  average  and  range  of  U values  for  the  stress-ratio  tests 

are  summarized  in  Figure  B-10(b).  It  is  apparent  that  closure  is  not  effective 

above  R = 0.15  for  these  conditions  and  not  particularly  dominant  at  R = 0. 

i 

Closure  effects  were  determined  at  K =10  and  20  ksi-in.  and  were  found  to 

max 

be  independent  of  K 

r max 


Conclusions 


The  experimental  data  generated  on  this  portion  of  the  research  program 
present  a refined  characterization  of  the  fatigue-crack-propagation  behavior 
of  T1-6A1-4V  plate  in  the  lower  range  of  AK  (i.e.,  10  to  20  ksi-in.^).  It  is 
apparent  that  the  behavior  is  influenced  by  processing  as  demonstrated  by  the 
differences  in  the  two  materials  evaluated.  The  response  to  overloads,  however, 
does  not  appear  to  be  greatly  different,  thereby  suggesting  that  the  retardation 
is  affected  by  the  plastic  zone  size  (or  related  residual  stresses)  rather  than 
the  processing  factors. 

In  testing  under  conditions  of  constant  stress-intensity-factor  range,  the 
crack-growth  curves  exhibit  far  more  linearity  (as  would  be  expected)  than  tests 
under  conditions  of  constant -amplitude  loading  and  confirm  the  basic  dominant 
role  of  AK  in  determining  rate  behavior.  However,  this  linearity  appears  to 
hold  only  for  crack  lengths  up  to  half  the  panel  width.  Divergence  from 
linearity  behond  this  point  indicates  that  either  the  finite-width  corrections 


are  inadequate  (which  is  unlikely)  or  that  other  mechanical  factors  come  into 
play. 

The  compliance  response  of  the  specimen,  as  well  as  the  crack-closure 
observations,  are  quite  consistent  with  current  concepts  of  crack  behavior.  This 
substantiates  the  validity  of  the  analytical  approach. 

Recommendations 


The  successful  completion  of  this  first  stage  of  analytical  and  experi- 
mental iteration  suggests  continuation  of  this  research  work.  It  is  recommended 
that  higher  ranges  of  AK  be  investigated  and  additional  stress  ratios  be  consid- 
ered. The  metallurgical  investigation  of  differences  in  the  two  materials  should 
also  be  pursued. 
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APPENDIX  C 


A SURVEY  OF  MICROSTRUCTURAL  INFLUENCES  ON  FATIGUE -CRACK- 
GROWTH  -RATE  DATA  WITH  SPECIAL  REFERENCE  TO  Ti-6A1-4V 


To  provide  a more  general  overview  of  the  mechanics-materials  interface 
relevant  to  fatigue-crack  growth,  this  program  was  augmented  with  a brief,  but 
general,  survey  of  relationships  observed  between  microstructural  features  and 
fatigue-crack  growth.  Although  special  reference  was  placed  on  Ti-6A1-4V,  other 
alloys  and  materials  were  included  to  provide  a more  definitive  interpretation 
of  microstructural  influences.  These  results  are  included  in  this  report  for 
general  information.  They  reflect  an  initial  basis  for  linking  microstructural 
effects  and  for  generalizing  the  evolving  model  to  other  materials.  It  should  be 
noted  that  the  International  System  of  Units  has  been  adopted  in  this  appendix 
to  be  consistent  with  the  majority  of  sources  which  were  included. 

Summary 

A review  of  the  published  literature  on  fatigue-crack  growth  suggests  that 
power-law  growth  in  Ti-6A1-4V  is  sensitive  to  microstructural  changes  which  result 
in  variations  in  the  fatigue  mechanism.  Microstructures  which  promote  secondary 
cracking  along  c*/8  interfaces  display  slow  growth  rates  while  microstructures 
which  promote  dimpled  rupture  display  fast  growth  rates.  Examples  of  similar 
effects  are  found  in  other  alloy  systems. 

Typically,  the  power-law  growth  exponent  in  this  alloy  is  3.4.  It  is  also 
suggested  that  the  power-law  regime  begins  at  AK  ~ 13  MNm"2,  coinciding  with  the 
lower  limit  of  striation  formation  on  the  fracture  surface.  The  upper  limit 
occurs  at  about  Kmax  = \ K^.  At  higher  growth  rates,  the  Forman  equation  appears 
to  be  adequate. 

The  normalized  stress-intensity  factor,  AK/E,  required  to  produce  a given 
growth  rate  in  Ti-6A1-4V  is  on  the  order  of  that  for  other  titanium-base  alloys, 
ferritic  steels,  martensitic  steels,  and  aluminum  alloys.  Austenitic  steels, 
which  deform  by  planar  slip  are  much  more  resistant  to  crack  growth  over  much  of 
the  stress-intensity  range  normally  encountered. 
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Background 


The  form  of  the  fatigue-crack-growth  curve  for  metals  and  alloys  is  quite 
well  established  and  is  typified  by  Figure  C-l  (Crooker  and  Krause  [1]*).  In  the 
curve  of  growth  rate  per  cycle,  Aa/An,  as  a function  of  stress-intensity  range, 

AK,  the  existence  of  three  separate  regions  is  evident.  The  central  power-law 
region  is  a well-defined  straight  line  on  a double- logarithmic  plot, 

ff  = c*(AK)m  . (C-l) 

The  major  emphasis  in  this  report  will  be  placed  on  the  region  described  by 
Equation  (C-l). 

The  individual  regions  of  Figure  C-l  roughly  correlate  with  differences  in 
fracture  appearance.  Yuen,  et  al  [2],  reported  the  quantitative  fractographic 
data  for  Ti-6A1-4V  plotted  in  Figure  C-2.**  The  moderate  values  of  AK,  corres- 
ponding roughly  with  the  region  of  power-law  crack-growth  behavior,  are  associ- 
ated with  a predominately  striated  fracture  surface.  These  striations  are  shallow 
lines  whose  spacing  is  close  to  Aa/An  for  Ti-6A1-4V  (as  well  as  for  a number  of 
other  alloys)  as  demonstrated  by  Bates  and  Clark  [3]  and  confirmed  by  Yuen,  et 
al  [2],  For  example,  Bates'  and  Clark's  data,  shown  in  Figure  C-3,  revealed 
that  the  striation  spacing  was  always  within  a factor  of  two  of  the  macroscopic 
growth  rate.  Figure  C-2  also  shows  that  at  both  the  high  and  low  ends  of  the 
AK  range,  other  microscopic  mechanisms  come  into  play.  The  high  values  of  AK 
(roughly  the  instability  region  of  Figure  C-l)  cause  the  introduction  of  the 
dimpled  rupture  mechanism,  characteristic  of  failure  under  monotonic  loading. 

At  the  other  end  of  the  scale,  the  low  values  of  AK  (roughly  the  threshold  region 
of  Figure  C-l)  are  associated  with  the  cleavage  mechanism. 

Power-Law  Fatigue-Crack  Growth 


Ti-6A1-4V 

Of  all  the  titanium-base  alloys,  Ti-6A1-4V  has  been  the  most  widely 
investigated.  The  results  in  Table  C-l  are  derived  from  data  of  10  different 


* References  are  listed  on  page  116. 

**  Throughout  this  appendix,  only  data  for  positive  stress  ratios  are  considered. 


Crack  Growth  Rate,  -r—  , /z.m/cycle 


FIGURE  C-l.  TYPICAL  FATIGUE-CRACK-GROWTH  CURVE  i 1] 


Striations  8 cleava 


ITY  RANGE  (2 


Best  Fit  to  Data 
Striation  Spacing  = growth  rate 
Striation  Spacing  = 1/2  growth  rate 
Striation  Spacing  = 2 growth  rate 


Striation  Spacing ,/xm 


FIGURE  C-3 . RELATION  BETWEEN  FATIGUE-CRACK-GROWTH  RATE  AND  SPACING  OF  STRIATIONS 
ON  THE  FRACTURE  SURFACE  OF  Ti-6A1-4V  L 3J 
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TABLE  C-l.  RESULTS  OF  POWER- LAW  ANALYSIS  OF  FATIGUE-CRACK 


GROWTH 

IN  Ti 

-6A1 

-4V 

Carve 

Number 

Yield 
Strength 
(MNnT3  ) 

Thick- 
ness 
(B,  mm) 

R 

M 

AK  (MNm*-) 

Normalized 
Actual  to  R = 0 ^ a ' 

Reference 

i 

875 

101.6 

0.0 

3.0 

24 

24 

Bates  & Clark  [3] 

2 

799(b) 

3.2 

0.33 

2.9 

20 

23 

Bucci,  et  al  [4] 

3 

944(b) 

0.2 

3.7 

22 

24 

4 

792^.b? 

_ _ 

0.0 

3.0 

31 

31 

Crooker  Si  Krause  [1] 

5 

320  b 

-- 

0.0 

3.3 

32 

32 

6 

903  b 

-- 

0.0 

3.3 

26 

26 

7 

985  b 

-- 

0.0 

3.3 

26 

26 

8 

944 

6.3 

0.1 

4.4 

18 

19 

Feddersen  Si  Hvler  [5] 

9 

944 

6.3 

0.4 

4.4 

16 

20 

10 

944 

6.3 

0.7 

4.2 

13 

19 

11 

834^) 

25.4 

0.08 

2.3 

22 

23 

Harrigan  [6] 

12 

841  b 

12.7 

0.30 

2.1 

17 

20 

13 

834  b 

25.4 

0.30 

1.6 

19 

22 

14 

1102(b) 

12.7 

0.30 

2.1 

14 

17 

15 

835  (b) 

10 

0.35 

3.1 

25(c) 

28 

Irving  Si  Beevers  [7] 

16 

863  b 

10 

0.35 

3. 1 

25  c 

28 

17 

1000  b 

10 

0. 12 

3. 1 

27  c 

28 

18 

1000  b 

10 

0.35 

3.1 

25  c 

28 

19 

1000 (b) 

10 

0.61 

3.1 

22  C 

27 

20 

935 

1.6 

0.  1 

2.5 

17 

18 

Jones  [8] 

21 

930 

12.6 

0.3 

3.9 

18 

21 

Paton,  et  al  [9] 

22 

937 

12.6 

0.3 

4.2 

22 

25 

23 

951 

12.6 

0.3 

5.4 

15 

18 

24 

1040 

12.6 

0.3 

5.4 

15 

18 

25 

974<d> 

0.8 

0.1 

3.8 

17 

18 

Pittinato  [10] 

26 

335  (b) 

0.33 

2.7 

2°(c) 

23 

Robinson,  et  al  [11] 

27 

950  b) 

-- 

0.33 

2.7 

20  c 

23 

28 

710 

6.4 

0.1 

3.5 

24 

23 

Thompson  [ 12  ] 

29 

772 

6.4 

0.1 

2.7 

28 

29 

30 

861 

6.4 

0.1 

3.5 

28 

29 

31 

903 

6.4 

0.1 

2.4 

20 

21 

32 

930 

10.2 

0.1 

5.2 

11 

12 

Thompson,  et  al  [13] 

33 

896 

10.2 

0.1 

4.4 

13 

14 

34 

10.2 

0.  1 

3.8 

15 

16 

35 

730!de! 

25.4 

0. 1 

3.1 

21 

21 

Yoder,  et  al  [14] 

36 

772  d) 

25.4 

0.1 

4.9 

24 

25 

37 

841 

25.4 

0.1 

4.2 

27 

28 

38 

916 

2.0 

0.0 

3.0 

22 

22 

Yuen,  et  al  [2] 

39 

916 

2.0 

0.1 

3.5 

21 

23 

40 

916 

2.0 

0.5 

3.2 

17 

22 

41 

916 

2.0 

0.7 

3.2 

15 

21 

42 

916 

2.0 

0.9 

3.2 

15 

23 

(a)  See  Figure  C-4. 

(b)  Variations  in  yield  strength  produced  by  variations  in  chemistry  and  processing. 

(c)  Extrapolated. 

(d)  Extra-low  intersitial  atom  content  (ELI). 

(e)  Ultra-low  interstitial  content,  0.06  wt.  percent  oxygen. 
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laboratories.  Only  room-temperature  tests  in  relatively  inert  environments 

(i.e.,  air,  inert  gas  atmospheres,  or  vacuum)  are  included.  The  data  are 

characterized  by  the  exponent  m of  Equation  ( C— 1 ) and  the  AK  value  required  to 

produce  a growth  rate  of  0.25  ^m/cycle  (10  x 10-3  in. /cycle);  the  latter  will 

be  denoted  AK7.  Three  of  the  investigations  involved  systematic  variations  in 

the  3tress  ratio,  R (R  = K . /K  ).  For  these  data,  AK  was  found  to  vary 

min  max 

linearly  with  R (see  Figure  C-4) , 

AK'  = ^ - 9R  , (C-2) 

where  AKg  is  the  value  of  AK7  for  R = 0.  Subsequent  comparisons  will  be  based 
on  AKq  derived  from  AK7  via  Equation  (C-2). 

Two  other  variables  are  considered  in  Table  C-l.  Processing  variations 
provided  a spread  of  yield  strengths  of  ~ 700  to  ~ 1100  MNm"‘ , while  thickness 
varied  by  two  orders  of  magnitude.  Thickness  appeared  to  have  no  significant 
effect  on  either  A Kg  or  m.  The  effect  of  yield  strength  is  not  clear,  because 
strength  changes  result  from  microstructural  changes  and  there  is  evidence  that 
fatigue-crack-retarding  microstructures  can  be  produced  (Barsom  [15]).  In 
particular,  secondary  cracks  have  been  observed  to  propagate  along  a/0  inter- 
faces inclined  to  the  propagation  direction  of  the  main  crack  in  beta-annealed 
Ti-6A1-4V  (Paton,  et  al  [9]).  In  turn,  these  interfaces  contain  a microstruc- 
tural phase  which  can  be  altered  by  heat  treatment  (Rhodes  and  Williams  [16]). 
It  is  possible  that  the  conflicting  reports  on  the  effect  of  oxygen  on  crack- 
growth  rate  (Thompson,  et  al  [13],  and  Yoder,  et  al  [14])  actually  reflect 
variations  in  the  grain  boundary  phase. 

Because  of  the  microstructural  effect,  it  is  not  possible  to  define  a 
"master"  fatigue-crack-growth  curve  for  Ti-6A1-4V.  However,  it  is  useful  to 
define  typical  behavior.  Accordingly,  average  values  taken  from  Table  C-l  are 
listed  in  Table  C-2,  along  with  results  for  the  other  alloy  classes  which  are 
discussed  below. 

A final  variable  to  consider  is  fracture  toughness.  The  only  systematic 
data  at  hand  are  a comparison  of  two  processing  treatments  resulting  in  similar 
yield  strengths  but  large  K£  differences  by  Crooker  and  Krause  [1]  which  result 
in  virtually  identical  power-law  growth  behavior. 


FIGURE  C-4.  EFFECT  OF  STRESS  RATIO  ON  STRESS  INTENSITY  TO  PRODUCE  A FATIGUE 
CRACK-GROWTH  RATE  OF  0.25  p,m/ CYCLE  IN  Ti-6A1-4V 


TABLE  C-2 


EVALUATION  OF  THE  PARAMETERS  DESCRIBING 
POWER-LAW  FATIGUE-CRACK  GROWTH 


(a)  Aa/An  (|j,m/cycle)  = 0.25  (■ 

for  Aa/An  = 0.25  y,m/cycle 
R = 0 via  Equation  (C-2). 

(b)  Average  values. 


■j  , where  AKq  is  the  stress-intensity  range 
Data  on  titanium-base  alloys  normalized  to 


Other  Titanium-Base  Alloys 


Data  for  several  other  titanium-base  alloys,  including  commercially  pure 
titanium,  are  given  in  Table  C-3.  This  allows  the  range  of  yield  strengths  to 
be  extended  to  200  - 1200  MNm_s  (Y/E  = 0.0011  - 0.006).  In  contrast,  these 
additional  alloys  appear  to  be  more  fatigue  resistant  than  Ti-6A1-4V  and  the  base 
metal  less  fatigue  resistant. 


Ferritic  and  Martensitic  Steels 


Several  significant  conclusions  have  resulted  from  investigations  on  steels 
The  data  in  Table  C-4  provide  results  from  selected  systematic  studies.  In  all 
cases,  particular  care  was  taken  by  the  authors  to  delineate  power-law  growth. 


Alloy 

Classification 

Elastic 
Modulus 
(MNnf  s ) 

(b) 

m 

AKo 

(MNtrf  T) 

A*o\ 

(10-4m2) 

Reference 

Titanium  Base 

Titanium  Metal 

110 

x 103 

3.4 

18 

1.6 

Table  C-3 

T1-6A1-4V 

110 

x 103 

3.4 

22 

2.0 

Table  C-l 

Other  Alloys 

110 

x 103 

3.3 

27 

2.5 

Table  C-3 

Steels 

Ferritic  and 

Martensitic 

200 

x 103 

2.65 

35 

1.75 

Table  C-4 

Austenitic 

200 

x 103 

6.1 

34  - 63 

1.70  - 3.15 

Thompson  ^ 1 

Aluminum  Base 

2024-T3 

68 

x 103 

2.6 

13  - 20 

1.9  - 2.9 

Hahn  & Simon 

7075-T6 

68 

x 103 

2.6 

10  - 15 

1.5  - 2.2 

If 

Others 

68 

x 103 

2.7 

12  - 16 

1.8  - 2.4 

fl 

(a)  Variations  in  yield  strength  produced  by  variations  in  grain  size 

(b)  Variations  in  yield  strength  produced  by  variations  in  processing 

(c)  Extrapolated. 

(d)  See  Figure  C-4. 

(e)  Estimated. 

(f)  Designated  "Plate". 


TABLE  C-3 . RESULTS  OF  POWER-LAW  ANALYSIS  OF  FATIGUE-CRACK  GROWTH 
FOR  TITANIUM  ALLOYS  OTHER  THAN  Ti-6A1-4V 


Yield 
Strength 
(MNciT3  ) 


Commercial  Purity 

V7(a) 
247  , v 

25° ^a) 

390(a) (c) (d) 

2 

0.35 

4.3 

16 

19 

Robinson  $•  Beevers  [19] 

2 

0.35 

4.3 

16 

19 

2 

0.35 

4.3 

16 

19 

Commercial  Purity 

220 

260 (a^ 
260  a 

.. 

0.07 

2.7 

17 

18 

Robinson,  et  al  [11] 

-- 

0.07 

2.7 

17 

18 

-- 

0.70 

2.7 

17 

23 

Ti  + 0.277.  02 

— 

6 

0.1 

3.1 

15 

16 

Thompson,  et  al  [13] 

Ti  + 0.57,  02 

-- 

6 

0.1 

3.1 

12 

13 

4A1-3MO-1V 

860SM 
860  b>e 
1000^’  l 
1000^’e) (f) 
1000^,e^ 
1000(b,e) 

1.25 

1.25 

0.05 

0.67 

2.5 

2.0 

28^ 

l7(c) 

33  (c> 

28 

23 

Boeing  [20] 

4.0 

0.05 

4.9 

33 

4.0 

0.67 

3.8 

19.  ^ 

33(c) 

25 

12.7 

0.05 

4.9 

33 

12.7 

0.67 

6.1 

19 

25 

6Al-6V-2Sn 

1095 fb) 
1126  b 

1038  b 
1194  J 
946  JJ 
1095  b 
U26  J 
1038 <b> 
1194 (b) 

(£) 

0.1 

2.6 

19 

20 

Amateau,  et  al  [21] 

(f) 

0.1 

2.6 

22 

23 

(f) 

0.1 

2.6 

25 

24 

(f) 

0.1 

2.6 

27 

28 

(f) 

0.1 

2.6 

30 

31 

(f) 

0.3 

2.9 

16 

19 

(f) 

0.3 

2.6 

21 

24 

(£) 

0.3 

2.6 

25 

28 

(f) 

0.3 

1.8 

28 

31 

7Al-2Nb-lTa 

758 

-- 

0.0 

5.0 

30 

30 

Crooker  [22] 

CN 

CO 

■ 

CN 

>% 

0) 

X 

rH 

3 

r*l 
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0 

CO 

C 

r^ 

tH 

rH 

mi 

X 

<D 

<V 

. 

0- 

— T 

U0 

uo 

X 

* 

rH 

CO 

CN 

CO 

>—i 

o 

03 

V 

E 

s 

Vi 

CQ 

U 

X 

O 

o 

«0 

iH 

V 

CO 

CO 

x 

X 

1—1 

u 

V4 

u 

u 

mi 

rH 

•H 

G 

CQ 

•H 

CQ 

-H 

CQ 

CQ 

CQ 

Q* 

CQ 

z 

<f  a\  in 

<o  co  co  co 


vC  O m co 
co  co  co  <o 


o ^ O o 

CO  CO  CO  CO  00 
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m m ^ co 
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X 

V 
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00 
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V 

CO 
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CO 

w 
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rH 

G 

G 
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Ml 
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M 
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c 
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a) 

00 
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03 
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O 

CO 
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rH 
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0 
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G 

G 

G 

G 

X 
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c 
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1 

1 

z 
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Li 

Li 

Li 

Li 

Li 

P" 

CQ 

CO 
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CN 
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3 
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Note  that  the  crack-growth  parameters  are  quite  reproducible  and  independent  of 
yield  strength  (Y/E  * 0.0006  - 0.007).  With  regard  to  the  mass  of  data  that 
do  not  display  values  close  to  those  of  Table  C-4,  Richards  and  Lindley  [19] 
and  Ritchie  and  Knott  [20]  point  out  that  those  data  tend  to  low  values  of  AKq 
and  high  values  of  m and  arise  when  the  micromechanism  departs  from  striation 
formation.  A similar  effect  can  be  noted  in  the  results  of  Paton,  et  al  [9], 
on  Ti-6A1-4V.  The  observed  high-growth  rates  reflect  the  intrusion  of  monotonic 
failure  modes  and  the  resultant  intrusion  of  the  instability  region  of  Figure  C-l 
and  the  accompanying  change  in  fracture  mode. 

A second  point  to  be  made  from  the  results  on  steel  is  that  the  effect  of 
R on  power-law  crack  growth  is  minimal  (Ritchie  and  Knott  [27]).  This  is  not 
the  case  for  titanium  alloys.  Thus,  it  does  not  currently  appear  to  be  possible 
to  generate  a universal  fatigue-crack-growth  law,  unless  separate  R variation 
is  defined  for  each  alloy  class.  Finally,  note  in  Table  C-2  that  AKg/E  is  about 
17  percent  higher  for  Ti-6A1-4V  than  for  steels,  suggesting  a slightly  higher 
inherent  crack-growth  resistance. 

Aluminum-Base  Alloys 

Hahn  and  Simon  [18]  have  reported  an  extensive  survey  of  the  literature 
on  high-strength  aluminum  alloys.  The  results  in  Table  C-2,  which  are  extracted 
from  their  survey,  are  for  dry  atmospheres  and  intermediate  frequencies  to  be 
most  comparable  to  the  data  of  other  systems.  These  results  appear  to  show  a 
metallurgical  effect,  with  AKq  being  somewhat  higher  for  2024-T3  than  for  7075- 
T6.  However,  Rosenfield,  et  al  [28],  have  shown  that  crack-growth  rates  in  7075 
can  be  reduced  to  those  characteristic  of  2024-T3  by  a proper  choice  of  heat 
treatment,  avoiding  the  toughness  minimum  associated  with  the  T-6  condition. 
Furthermore,  Garrett  and  Knott  [29]  have  evidence  that  intrusion  of  monotonic 
fracture  modes  increases  both  the  m value  and  the  apparent  effect  of  R.  Thus, 
it  is  possible  that  the  "real"  behavior  of  aluminum  alloys  is  closer  to  that  of 
steel  than  generally  believed. 

Austenitic  Steels 


The  only  major  difference  in  power- law  behavior  among  the  alloys  surveyed 
was  found  in  austenitic  steels.  Over  most  of  the  range,  growth  rates  for  these 
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materials  are  significantly  slower  than  for  all  of  the  other  alloys.  The 
results  cited  in  Table  C-2  are  from  Thompson's  [17]  survey.  He  points  out  that 
the  large  exponent  is  observed  despite  the  presence  of  striations.  Bathias  and 
Pelloux  [25]  report  that  the  curve  of  stristion  spacing  versus  AK  is  the  same 
for  austenitic  and  maraging  steels.  However,  the  striations  in  maraging  steels 
are  aligned  normal  to  the  crack  front  while  those  in  austenitic  are  at  various 
angles  in  different  grains.  They  suggest  that  the  crack  front  advances  piece- 
meal along  its  length,  accounting  for  the  discrepancy  in  macroscopic  and  micro- 
scopic growth  rates. 

These  observations  have  implications  for  development  of  models  of  crack 
growth.  Austenitic  steels  are  an  extreme  example  of  planar  slip  materials-- 
dislocation  lines  tend  to  remain  straight  and  not  to  cross  slip.  In  contrast, 
the  other  materials  surveyed  cross  slip  extensively  and  form  dislocation  tangles. 
This  well-known  difference  in  monotonic  deformation  patterns  is  also  observed 
after  cyclic  straining  (Abdel-Raouf , et  al  [30]).  The  extreme  change  in  micro- 
scopic deformation  behavior  apparently  does  not  affect  the  microscopic  growth 
rate  (i.e.,  striation  configuration)  leading  to  the  probability  that  highly 
artificial  dislocation  models  of  crack-tip  deformation  patterns  may  be  suffi- 
cient to  analyze  crack  growth. 

Upper  and  Lower  Limits  of  the  Power-Law  Range  in  Ti-6A1-4V 
Threshold  Region 

As  yet,  there  is  not  sufficient  evidence  to  define  threshold  behavior.  On 
the  particularly  important  question  of  the  effect  of  strength  and  microstructure, 
the  evidence  is  conflicting.  Irving  and  Beevers  [7]  find  that  microstructural 
changes  strongly  affect  Aa/An  in  the  threshold  region,  while  Crooker  and  Krause 
[1]  come  to  the  opposite  conclusion.  The  situation  is  compounded  because  this 
is  the  region  of  very  low  growth  rates  where  experimental  errors  are  expected 
to  be  greatest.  As  a result,  it  is  not  yet  possible  to  define  a general  fatigue- 
crack-growth  law  for  rates  below  the  power-law  region  for  Ti-6A1-4V.  However, 
Weiss  and  Lai  [31]  have  suggested  a three-parameter  curve  based  on  Neuber 
analysis  which  appears  to  be  a reasonable  basis  for  further  development. 
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V 


The  data  in  Table  C-5  suggest  that  the  transition  from  threshold  behavior 

to  power-law  growth  occurs  in  the  range  0.01  - 0.4  jj.ni/cycle,  corresponding  to 

AK  values  of  ~ 8 - 20  MNnf^.  This  corresponds  to  the  cleavage-striation 

mechanism  change  observed  by  Yuen,  et  al  [2].  It  is  suggested  that  an  inter- 

2 

mediate  value  of  AK  = 13  MNnf  2 be  tentatively  adopted  as  the  lower  limit  of 
power- law  growth. 

Instability  Region 


As  with  threshold  behavior,  it  is  not  yet  possible  to  make  any  definite 
conclusions  regarding  the  onset  of  instability.  The  appearance  of  dimpled  rup- 
ture on  the  fracture  surface  suggests  a connection  with  Kc  and  there  is  an 
indication  that  the  AK  values  associated  with  the  upper  cutoff  of  the  power-law 
region  are  higher  in  tougher  alloys.  To  account  for  crack  growth  in  the 
instability  region,  it  is  suggested  that  the  empirical  equation  of  Forman, 
et  al  [33]  be  used  since  it  has  been  applied  successfully  to  Ti-6A1-4V  (Fedder- 
sen  and  Hyler  [5]).  A modified  form  of  Equation  (C-l)  is 


A£  = 

An 


aAK 


(1_R)  AK  - 1 


(instability  region) 


(C-3) 


The  upper  cutoff  of  the  power-law  region  is  then  found  by  equating  Equations 
(C-l)  and  (C-3)  so  that  the  transition  condition  becomes 


AK  1 - R 

K 2 
c 


(C-4) 


The  available  data  in  Table  C-5  are  in  fair  agreement  with  Equation  (C-4)  so 
that  it  can  be  considered  a reasonable,  if  not  exact,  criterion. 


Discussion 


Research  in  fatigue-crack  growth  over  the  past  few  years  has  revealed  that 
microstructural  variations  can  play  a significant  role.  Depending  on  the  AK 
level  and  on  the  composition/processing  history,  at  least  four  different 
mechanisms  are  operative  in  Ti-6A1-4V:  cleavage,  dimpled  rupture,  striation 

formation,  and  secondary  cracking.  While  not  enough  information  is  available 


on  cleavage,  it  appears  that  the  growth  rates  can  be  improved  by  promoting 
secondary  cracking  and  inhibiting  dimpled  rupture.  Also,  it  is  reasonable  to 
conclude  that  any  crack-growth  equation  must  take  into  account  the  micromecha- 
nism. For  example,  the  value  of  K should  enter  into  the  equation  in  a differ- 

nicix 

ent  way  when  the  fracture  surface  displays  dimpled  rupture  than  when  it  is 
completely  striated. 

Because  of  these  microstructural  variations,  it  is  not  yet  possible  to 
define  a "master  curve"  for  fatigue-crack  growth  in  Ti-6A1-4V.  Instead,  equa- 
tions of  "typical  behavior"  are  offered  in  Table  C-6.  The  equations  are  useful 
for  making  preliminary  estimates  of  the  behavior  of  off-the-shelf  pieces  and 
provide  a standard  to  decide  whether  a given  microstructure  displays  outstand- 
ingly good  or  outstandingly  poor  behavior.  Furthermore,  it  should  be  emphasized 
that  the  microstructural  effect  is  sufficiently  large  that  as  of  now  actual 
crack-growth  measurements  are  the  only  satisfactory  method  of  accurately  deter- 
mining growth  rate  data  in  a given  piece  of  this  alloy. 

Despite  this,  the  relative  insensitivity  of  growth  rates  via  the  striation 
mechanism  to  normalized  stress  intensity,  AK0/E  (Table  C-2),  remains.  This 
well-known  phenomenon  provides  an  unsolved  problem  and  suggests  some  difficulties 
in  attaining  major  improvements  in  growth-rate  resistance. 
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